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Abstract 

I review recent progress on understanding QCD dynamics involved in exclusive B meson decays. 
Different frameworks, including light-cone sum rules, QCD factorization, perturbative QCD, soft- 
collinear effective theory, light-front QCD, are discussed. Results from lattice QCD are quoted for 
comparison. I point out the important issues in the above QCD methods, which require further 



(N 



investigation. 



2 1 Introduction 

> 



We are now in the era of B physics. B fatories at KEK and SLAC have collected about 80 data, 
\Q | based on which we are not only able to probe the origin of CP violation, but to explore rich QCD 
7—1 ; dynamics involved in exclusive B meson decays. As announced in |T|, the Kobayashi-Maskawa (KM) 
^ ■ ansatz 2 J for CP violation is more or less certain with the consistent measurements of sin 20! (or sin 2/3) 
O ■ from Belle and BaBar, 0! being one of the unitarity angles. The results are also in agreement with other 
^ ! indirect determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. In this article I 
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will forcus on the latter subject. It will be realized that exclusive B meson decays provide a unique 



Q_i' field, in which QCD theories with controllable theoretical uncertainty can be developed. It turns out 



o ■ that these theories are simpler than those for charm and kaon physics. This field has attracted wide 
■ attention, and tremendous progress has been made recently. 

Within the KM ansatz, the source of CP violation is organized in the form of a unitarity triangle. 
.£h ! On one hand, we overconstrain the unitarity triangle as much as possible, and on the other hand, look 
^ \ for possible discrepancies, which could reveal signals of new physics beyond the Standard Model. The 
& ' angle 0i can be extracted from the CP asymmetry in the B — > J/ipKg decays [3], which arises from the 
B-B mixing. Through similar mechanism, the B® — ► 7r + 7r _ decays are appropriate for the extraction 
of the angle 02 (or a). However, these modes contain both tree and penguin contributions, such that 
the extraction suffers uncertainty. Strategies have been proposed to handle this penguin pollution, the 
best of which is known to be the isospin analysis [3]. Unfortunately, this strategy is difficult in practice, 
because of the small B® — > 7r°7r branching ratios. The angle 03 (or 7) can be determined from the 
decays B — > Kit jSJ |H1 13 El , which are obviously also plagued by the similar penguin-tree interference. 

We can move forward by learning how to estimate hadronic matrix elements involved in exclusive 
B meson decays. For this purpose, symmetries of strong interaction have been postulated to relate 
amplitudes among different modes. For example, the penguin-over-tree ratio \P/T\ helps the extraction 
of 02 from the CP- violating observables in B® — ► tt + it~ <S77(3) flavor symmetry and plausible 
dynamical assumptions were then employed to fix \P\ through the CP-averaged B ± — > Kn^ 1 branching 
ratio ^0]. The information of |T| can be obtained from the B — > nlu data. Another strategy is to 
apply the [/-spin flavor symmetry to the B% — > 7r + 7r~ and B s — ► K + K~ modes jllj . from which the 
penguin amplitudes are determined. However, the above symmetries are in fact not exact, and it is 
not clear how to control theoretical uncertainties from symmetry breaking. As an alternative, one 
searches for the special modes, in which relations among decay amplitudes allow the elimination of 
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A =C {L' U , L> , Uy } amplitudes \12\ I13|. which receive only tree contributions, being the Cr-even 
eigenstate of the neutral .D-meson system. However, this strategy, due to the squashed triangles, is 
experimentally difficult [H]. The modes B d -> K*°{D°, D°, D° + } [IS] and B± -> Df{D ,D ,D° + } jTH], 
providing more equilateral triangles, then deserve further feasibility studies. 

The above discussion indicates that it is necessary to have deeper understanding of QCD dynamics 
in exclusive B meson decays and control of hadronic uncertainties |17j . The b quark mass mj, much 
larger than the QCD scale Aqcd, renders such an attempt possible: relevant hadronic matrix elements 
can be evaluated as an expansion in the strong coupling constant a s {mb) and in the ratio Aqcd/^6- The 
approaches based on this heavy quark expansion include light-cone QCD sum rules (LCSR) [TS1 |2U] . 
light-front QCD (LFQCD) [211 122, QCD factorization (QCDF) [23], and perturbative QCD (PQCD) 
[2U [213 12E] • Soft-collinear effective theory (SCET) provides a systematic framework, in which the above 
expansion can be constructed in a simple and formal way [27JI2E]- Lattice QCD is complementary to 
the above methods, whose results will be quoted for comparison. In this article I will explain the basic 
ideas behind the various QCD theories, and review their applications to typical, such as semileptonic, 
radiative and nonleptonic, exclusive modes. That is, I emphasize the methodology, instead of the survey 
of all decay channels. 

To be specific, I will not discuss the strategies to constrain the CKM matrix elements from exper- 
imental data. For nice reviews of this topic, refer to |2H1 HUB EH] • I will not explore analyses relying 
on symmetries of strong interaction, such as the SU(3) flavor symmetry. Recent works along this vein, 
which have taken into account symmetry breaking effects, can be found in [3*2*1 [331 EH] • The status of 
the important CKM global fitting have been presented in [!""""| ESI- To demonstrate the applications 
of the QCD methods, I will consider only B u & meson decays as an example. The subjects related to 
the B s and B c mesons and to heavy baryons, including their polarization effect, will be dropped. The 
perspectives for investigating B c mesons at LHCb have been surveyed in [37]. For a similar reason, I 
skip the applications to decays into baryons and into tensor mesons j3H[. Studies of three-body B meson 
decays are still at the early stage [3~""1 """U] . and will be reserved for a future review. I will not touch 
supersymmetric topics in B physics either, which are too much for this article. For a recent relevant 
review, refer to (HI H2] ■ 

In Sec. 2 and Sec. 3 I briefly explain two types of factorization theorems, collinear and kx factoriza- 
tions, which are the fundamental concepts of most of the QCD theories. The ideas and results derived 
from the various QCD methods are reviewed in Sec. 4 for semileptonic and radiative decays, and in 
Sec. 5 for two-body nonleptonic decays. Charmed decays are discussed in Sec. 6. Other miscellaneous 
topics are collected in Sec. 7. Section 8 is the conclusion. 

2 Collinear Factorization 

Most of QCD methods rely on some sorts of factorization theorems. For example, QCDF is a general- 
ization of collinear factorization theorem to exclusive B meson decays. In LCSR collinear factorization 
applies to final-state hadron bound states, which are then expanded in terms of parton Fock states 
characterized by different twists, kx factorization theorem is the basis of the PQCD approach, which 
is more appropriate in the end-point region of parton momentum fractions. SCET for the kinematic 
region with energetic final-state hadrons, is equivalent to collinear factorization theorem, but operated 
at the operator and Lagrange level. I will compare the factorization of high-energy QCD processes 
derived from perturbation theory and from SCET. I then discuss the double logarithmic resummation, 
which is required for applying collinear factorization theorem to semileptonic B meson decays. 

2.1 Factorization Theorem 

I first review collinear factorization theorem for exclusive processes developed around 80's [""3~l [33J ["""3 
HHIEIIIH]- In this theorem nonperturbative dynamics of a high-energy QCD process, characterized by 
a large scale Q, either cancel or is absorbed into hadron distribution amplitudes. The remaining part, 
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Figure 1: Lowest-order diagrams for tc-/* — > 7 (B — > 7/^), where the symbol x represents the virtual 
photon (weak decay) vertex. 



being infrared finite, is calculable in perturbation theory. A physical quantity is then expressed as the 
convolution of a hard-scattering kernel with the distribution amplitudes solely in parton momentum 
fractions. A distribution amplitude, though not calculable, is universal, i.e., process- independent. With 
this universality, a distribution amplitude determined by nonperturbative means, such as QCD sum rules 
and lattice QCD, or extracted from experimental data, can be employed to make predictions for other 
processes involving the same hadron. Contributions of different orders in a s and powers in 1/Q can be 
included systematically. 

Nonperturbative dynamics is reflected by infrared divergences in radiative corrections, whose fac- 
torization leads to distribution amplitudes at the parton level. Factorization of the above infrared 
divergences needs to be performed in momentum, spin, and color spaces. Factorization in momentum 
space means that a hard kernel depends on the loop momentum of a soft or collinear gluon, which 
has been absorbed into a distribution amplitude, only through the parton momentum fraction. Fac- 
torization in spin and color spaces means that there are separate fermion and color flows between a 
hard kernel and a distribution amplitude, respectively. I take the simple process rcy* — > 7 as an ex- 
ample to demonstrate the proof of collinear factorization thheorem based on perturbation theory. The 
collinear factorization of this process has been proved in [13], but in the axial (light-cone) gauge A + = 0. 
In this gauge factorization automatically holds and the analysis is straightforward, because collinear 
divergences exist only in two-parton reducible diagrams. The pion distribution amplitude has been 
constructed from 7*7 — > it in the framework of covariant operator product expansion |H| The 
factorization of tcj* ir has been proved in jlH] based on the Zimmermann's "A-forest" prescription 
49J, which involves complicated diagram subtractions. 

Below I will adopt a simple proof proposed in [20]. To achieve factorization in momentum, spin, 
and color spaces, one needs the eikonal approximation for loop integrals in leading infrared regions, the 
insertion of the Fierz identity to separate fermion flows, and the Ward identity to sum up diagrams 
with different color structures. Under the eikonal approximation, a soft or collinear gluon is detached 
from the lines in a hard kernel and in other distribution amplitudes. The Fierz identity decomposes 
the full amplitude into contributions characterized by different twists. The Ward identity is essential 
for proving factorization theorem in a nonabelian gauge theory. The soft divergences exist in exclusive 
B meson decays, which should be factorized into a B meson distribution amplitude. The derivation 
in [HII] is explicitly gauge-invariant, and appropriate for both the factorizations of soft and collinear 
divergences, compared to those in the literature. 

The momentum Pi of the pion and the momentum P 2 of the outgoing on-shell photon are chosen, 
in light-cone coordinates, as 



Let e denote the polarization vector of the outgoing photon, which contains only the transverse com- 
ponents. Consider the kinematic region with large Q 2 = —q 2 , q = P 2 — P\ being the virtual photon 
momentum, where perturbative expansion is reliable. The lowest-order diagrams are displayed in Fig. ^ 
Assume that the on-shell quark and antiquark carry the fractional momenta xP\ and xPi, respectively, 
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Figure 2: 0(a s ) radiative corrections to Fig. (TJ^a) 



with x = 1 — x. Figure HJa) gives the parton-level amplitude, 

QW(x) = -ie 2 q(xP 1 ) 4 %Zl£)2 W(* P i) 



(2) 



The analysis for Fig. ^b) is the same. The internal quarks are regarded as being hard, i.e., being 
off-shell by 0(Q 2 ). 

The factorization in the fermion flow is achieved by inserting the Fierz identity, 

hjhk = -rhkhj + -r(l5)ik(l5)ij + ](ja)ik(j a )ij + -;{lsl a )ik{lals)ii + l(v aP )ik((T*p)ij , (3) 



4 



where I represents the identity matrix, and a a p is defined as a a p = i[7 a ,7^]/2. Different terms in the 
above identity lead to contributions of different twists. Equation (J2J) is then factorized into 



(4) 



where the functions, 



2 tr(4 P 2 7^ Alb) 



APf 



g(xPi)7 5 fiq{xPi) 



2xP 1 ■ P 2 



(5) 



with the dimensionless vector n = (0, 1, Ot) on the light cone, define the lowest-order distribution 
amplitude and hard kernel in perturbation theory, respectively. For the momenta chosen in Eq. (J!}, 
only the pseudo-vector structure 70,75 with a = + survives, as it is contracted with the hard kernel 
to form the factor /Pi 75 in Eq. (jSJ). This piece of contribution is of leading twist (twist 2). For other 
processes, such as the pion form factor, higher-twist structures survive, but the analysis is the same 

EH- 

There are two types of infrared divergences in radiative corrections, soft and collinear. In the soft 
region and in the collinear region associated with the pion momentum P\, the components of a loop 
momentum I behave like 



P = (/+, T, 1 T ) ~ Q(X, A, A) , F ~ Q(l, A 2 , A) 



(6) 
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diminishes as A , and the corresponding loop integrand may diverge as 1/A . As the phase space tor 
loop integration vanishes like dH ~ A 4 , logarithmic divergences are generated. 

We identify the infrared divergences in the 0(a s ) corrections [HH E21 EB] to Fig. da), which are 
shown in Fig. |21 Self-energy correction to the internal quark, giving a next-to-leading-order hard kernel, 
is not included. The collinear factorization formula is written as the convolution over the momentum 
fraction £, 

g\ x \x) = Jd^\x,m {0) (O + <P {0) (x)n?\x), (7) 

The above expression, with the 0(a s ) distribution amplitudes <f>i(x, £) specified, defines the 0(a s ) hard 
kernels ri^\x)^ which do not contain collinear divergences. This procedure is referred to as matching 
the effective theory to the full thoery in the determination of Wislon coefficients in SCET. It is now 
obvious why an arbitrary x is considered for the parton- level diagrams in Figs. d an d 121 one can obtain 
the functional form of TQ (x) in x. 

Figures E[a)-|2tc) are the two-particle reducible diagrams with the additional gluon attaching the two 
valence quarks of the pion. It has been known that soft divergences cancel among these diagrams [50J. 
The reason for this cancellation is that soft gluons, being huge in space-time, do not resolve the color 
structure of the pion. Collinear divergences in Figs. I21^a)-I2tc) do not cancel, since the loop momentum 
flows into the internal quark line in Fig. Efb), but not in Figs. Efa) and|2fc). To absorb the collinear 
divergences, one introduces a nonperturbative distribution amplitude. The factorization of the above 
diagrams is achieved by inserting the Fierz identity. For example, one obtains, from Fig. Efb), the 
0(a s ) distribution amplitude, 

<t>b {x,0 = —p-r- / 77m^( xP i) T~F, IW /-n I 7x9,9 q(xP 1 )8[£-x + 



4P X + J (2vr) 4 ^ 17 {xP 1 -l) 2 {xP 1 + l) 2 l 2 * v 17 V Pf ) ' v ' 

with Cf = 4/3 being a color factor. <pjp contains the collinear divergence, because the integrand in 
Eq. (jHJ) diverges as 1/A 4 . The dependences on l~ and on It in 7i^°\ being subleading according to 
Eq. ©, have been neglected. 

In the collinear region of Fig. I2fd), the following approximation for part of the loop integrand holds, 

(fa-x frWifr-x ft+fi « 2P» / P 2 , (9) 

where the l~ and It terms, being power-suppressed compared to P 2 j have been dropped. The factor- 
ization of the collinear divergence from Figs. Efd) requires the further approximation for the product 
of the two internal quark propagators [50J, 



IPX n u r 1 1 



(P 2 -xP l ) 2 (P 2 ~xP 1 + l) 2 n-l 



(P 2 - xPi) 2 (P 2 - xPx + If 



(10) 



which is an example of the Ward identity. Similarly, the power-suppressed terms, such as I 2 and xP\ ■ I, 
have been neglected . The numerator 2P% comes from Eq. (|9*jl. and the factor n u /n ■ I is exactly 
the Feynman rule associated with a Wilson line. Therefore, the appearence of the Wilson line is a 
consequence of the Ward identity. 

The first (second) term on the right-hand side of Eq. (|1U|) corresponds to the case without (with) the 
loop momentum I flowing through the hard kernel. Hence, the extracted 0(a s ) distribution amplitude 
is written as 



,(i), c s -ig A C F f dH 5 x /Px+ / ,_ p ,l n" 





I 2 n-l 

(11) 
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side ol Eq. (jlOj) . ihe factorization of the distribution amplitude from r igs.|2Je) is performed m a similar 
way. 

The above 0(a s ) distribution amplitudes can be reproduced by the 0(a s ) terms of the following 
nonlocal matrix element with the structure 75 /i sandwiched, 



/dy + — py 

— — e~ l ^ PlV (0|g(?/ _ )75 fiP exp — ig / dzn ■ A(zn) 
Z7T JO 



qiOMxPjqim)) ■ (12) 



The notation V means the path ordering for the Wilson line. The integral over z contains two pieces: 
to generate the first term in Eq. (fTTj) . z runs from to 00; to generate the second term, z runs from 
00 back to y~ . The light-cone coordinate y~ 7^ corresponds to the fact that the collinear divergences 
in Fig. El do not cancel. The Wilson line along the light cone collects collinear gluons in irreducible 
diagrams. By expanding the quark field q{y~) and the Wilson line into powers oiy~, the above matrix 
element can be expressed as a series of covariant derivatives (D + ) n q(0), D = d — igA, implying that 
Eq. (fT^j) is gauge- invariant. 

I then review the all-order proof of leading-twist collinear factorization theorem for the process 
7T7* — > 7, and justify the definition of the parton- level distribution amplitude in Eq. (fT2*|) . The proof 
is performed in the covariant gauge, in which collinear divergences also exist in two-particle irreducible 
diagrams. It has been mentioned that factorization of a QCD process in momentum, spin, and color 
spaces requires summation of many diagrams, especially at higher orders. Hence, the diagram summa- 
tion must be handled in an elegant way. For this purpose, one employs the Ward identity, 

l ll G' i (l,k 1 ,k2,--',k n )=0, (13) 

where represents a physical amplitude with an external gluon carrying the momentum I and with n 
external quarks carrying the momenta k\, k 2 , ■ ■ •, k n . All these external particles are on mass shell. The 
Ward identity can be easily derived by means of the Becchi-Rouet-Stora (BRS) transformation [5"l*j . 

Factorization theorem can be proved by induction. The factorization of the 0(a s ) collinear diver- 
gences associated with the pion has been worked out in Eq. (J7|). Assume that factorization theorem 
holds up to 0(af), 

gV)( x ) = I d^\x,on^\0 , j = 1, ■ ■ ■ , N , (14) 

i=0 J 

where </>W(a;,£) is given by the 0{a\) terms in the perturbative expansion of Eq. (H2J), and H^H) 
stands for the 0{a{~ % ) infrared-finite hard kernel. It will be shown that the 0(a^ +1 ) diagrams Q( N+1 ) 
is written as the convolution of the O(af) diagrams Q^ N > with the 0(a s ) distribution amplitude by 
employing the Ward identity in Eq. (|T3J). 

Look for the gluon in a complete set of 0(a^ +1 ) diagrams Q^ N+l \ one of whose ends attaches the 
outer most vertex on the upper quark line in the pion. Let a denote the outer most vertex, and j3 
denote the attachments of the other end of the identified gluon inside the rest of the diagrams. There 
are two types of collinear configurations associated with this gluon, depending on whether the vertex 
(3 is located on an internal line with a momentum along Pi. The quark spinor adjacent to the vertex 
a is q(xPi). If f3 is not located on a collinear line along Pi, the component 7 + in 7° and the minus 
component of the vertex (3 give the leading contribution. If j3 is located on a collinear line along Pi, (3 
can not be minus, and both a and (3 label the transverse components. This configuration is the same 
as of the self-energy correction to an on-shell particle. 

According to the above classification, one decomposes the tensor g a/3 appearing in the propagator 
of the identified gluon as 

XX ( , X X ^\ flK\ 

9a(i = ——j- - OarOpT + \9aP ~ — J + a T0f3T I • (15) 

The first term on the right-hand side extracts the first type of collinear enhancements, since the light- 
like vector n a selects the plus component of 7°, and the dominant component ln=- in the collinear 
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Figure 3: (a) Ward identity, (b) Factorization of Q^ N+l \ 



region selects the minus component of the vertex (3. The components lp=+,T do not change the collinear 
structure, since they are negligible in the numerators compared to the leading terms proportional to 
Pi or P 2 ~. This can be confirmed by contracting lp to Fig. EJd), from which Eq. (JTUj) is obtained. The 
second term extracts the second type of collinear enhancements. The last term does not contribute a 
collinear enhancement due to the equation of motion for the valence quark. We shall concentrate on 
the factorization of Q^ N+1 ^ corresponding to the first term on the right-hand side of Eq. (JT5j) . and the 
factorization associated with the second term can be included following the procedure in |5D] . 

The contraction of lp hints the application of the Ward identity in Eq. (fT3*j) to the case with two 
external on-shell quarks. Those diagrams with Figs. Ufa) and 12(b) as the 0(a s ) subdiagrams are 
excluded from the set of G^ N+1 ^ as discussing the first type of collinear configurations, since the identified 
gluon does not attach a line parallel to P\. Consider the physical amplitude, in which the two on-shell 
quarks and the on-shell gluon carry the momenta £Pi, xP\ and I, respectively. Figure a), describing 
the Ward identity, contains a complete set of contractions of lp, since the second and third diagrams 
have been added back. The second and third diagrams in Fig. Of a) lead to 

p / A(W = 7 p 1 , (/- e Pi + 1 AM£A) = -qm) , 

IpqtxPr)^—^ = -q{xPi) , (16) 

x /Pi- / 

respectively. The terms q(£Pi) and q(xPi) at the ends of the above expressions correspond to the 0(a^) 
diagrams. 

Figure 0(b) shows that the diagrams Q^ N+1 ^ associated with the first term in Eq. (|15J) are factorized 

into the convolution of the parton-level O(a^) diagrams with the 0(a s ) collinear piece extracted 
from Fig. EKd). The factor n a /n ■ I from the collinear replacement in Eq. (|15J) is exactly the Feynman 
rule associated with the Wilson line in the direction of n, represented by the double line. The first 
diagram means that the gluon momentum does not flow into G^ N \ while in the second diagram the 
gluon momentum does. The similar reasoning applies to the identified gluon, one of whose ends attaches 
the outer most vertex of the lower antiquark line. Substituting Eq. (fUjl into on the right-hand 

side of Fig. HHb), and following the procedure in [SU], one arrives at 

N+l 

<3 {N+l \x) = E / d^\^0n {N+1 -H0 , (i7) 

i=0 J 

with the infrared- finite 0(a^ +1 ) hard kernel l~i^ N+1 \ Equation (|17j) implies that all the collinear en- 
hancements in the process rc-f* — > 7 can be factorized into the distribution amplitude in Eq. (J 12)) order 



7 



Equation ()12j) plays the role oi an mirared regulator tor the parton-ievel diagrams. A hard kernel 
H(x,Q 2 , fi) can then be regarded as the regularized parton-level diagrams, where the dependences on 
the momentum transfer Q and on the factorization scale \x have been made explicit. Different choices of 
fi correspond to different factorization schemes. Since the parton-level diagrams with on-shell external 
particles and Eq. ()12|) are both gauge- invariant, the hard kernel is gauge-invariant. After determining 
the gauge-invariant infrared-finite hard kernel, one convolutes it with the physical two-parton pion 
distribution amplitude, whose all-order gauge-invariant definition is given by 



^(x,ii) ^ij^-e ixP ? y (0\q(y ) 75 fiVexp -ig jT <lin--A(zn) 



Q 



7T 



(Pi)) ■ 



The valence-quark state \q(xPi)q(xPi)) has been replaced by the pion state \tt(Pi)), and the pion decay 
constant /„. has been omitted. Equation ([18)1 can also be derived in SCET as argued in the next 
subsection. The 717* — > 7 scattering amplitude is then expressed as the convolution over the parton 
momentum fraction x, 



M(Q 2 



dx<j) n (x, fj)TC(x, Q 2 , jj) 



(19) 



Hence, predictions derived from collinear factorization theorem are gauge-invariant and infrared- finite. 



2.2 Soft- Collinear Effective Theory 

Final-state hadrons in exclusive B meson decays may carry energy E of O(ms), tub being the B meson 
mass, which is much larger than Aqcd- These processes can be analyzed in the collinear factorization 
framework discussed in the previous subsection. To study the collinear factorization at the operator 
and Lagrange level, SCET has been developed j23 EHl IEH1 EE] • After integrating out short-distance 
fluctuations characterized by the invariant mass p 2 ^> (EX) 2 , which appear in Wilson coefficients, 
long-distance fluctuations are then described by new effective degrees of freedom. SCET then exhibits 
symmetries in the large energy limit, such as the reduction of spin structures, helicity constraints, and 
collinear gauge invariance, which apply to the new effective fields. Power corrections in SCET are 



included in terms of the small parameter A = Aqcd/E (or A = ^JAqcd/E) [57| EE1 EH| • For a recent 
review, refer to [60J. 

The effective fields contain collinear quarks and gluons (£ n , p , A£ ), massless soft quarks and gluons 
(q s , A%), and massless ultrasoft (usoft) quarks and gluons (q us , A£ s ). The collinear fields, labelled by 
the light cone direction n = (1,0, Ot) and their momentum p, come from the phase redefinitions, 

^(x) = ^e- !p -> n ». (20) 
p 

Derivatives on the collinear fields, 9 /i n)P (x) ~ (EX 2 )cf) nt p(x), pick up only the small scale. The large 
momenta are picked up by introducing the label operator, V£ n<p = (n-p)£, n ,p- Similarly, the operators 
Vt and P M are defined to pick up the 0(A) labels of the collinear and soft fields, respectively. In the 
discussion in Sec. 12.11 the usoft fields can be regarded as the leftover pieces with the collinear or soft 
dynamics being factorized out. That is, the usoft fields, due to their slow variation, appear as the 
background fields to the collinear or soft ones. 

The collinear Wilson line W[n ■ A n ^ q ] and soft Wilson line S n [n ■ A s ] are induced to preserve gauge 
invariance [^HUSHl- For the former, the explicit expression is given by 

W = exp[-ln-A n , q (x)} , (21) 

perms ' 

which is equivalent to that in Eq. (J18)) derived from perturbation theory. The meaning of a collinear 
Wilson line in the definition of the pion distribution amplitude has been emphasized in Sec. 12.11 The 
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collmear 




« P , ^.p) ~ (1,A 2 ,A) 


A', VT ~ A u 
P£ ~ A 


soft 


p"~(A,A,A) 


~ A a / 2 
< P ~ A 


V ~ A 


usoft 


fc^ ~ (A", A 2 , A 2 ) 


5w ~ A 3 
~ A 2 

us 





Table 1: Power counting rules for momenta, fields, momentum label operators (V,Vj>, V 11 ), and 
collinear, soft and usoft Wilson lines (W, S n , Y n ) in SCET. 

usoft Wilson line Y n [n ■ A us ] is introduced by the further field redefinitions = Y n $ffl p and A n>q = 

Assuming the action for the kinetic terms in SCET to be of 0(A°), the scaling of each effective field 
in A can be defined straightforwardly. The power counting rules for the momenta, fields, momentum 
label operators, and collinear, soft and usoft Wilson lines are summarized in Tabled 071 128j. It is 
found that the scaling of the collinear and soft momenta in SCET is the same as that defined in Eq. (JH)l. 
The usoft momentum in the framework based on perturbation theory will appear as discussing the 
factorization of soft divergence from exclusive B meson decays in Sec. 13.21 

The leading-order Lagrangians for (u)soft light quarks and gluons are the same as in QCD. For heavy 
quarks h v labelled by the velocity v, we have the heavy quark effective theory (HQET) Lagrangian |5Tj . 

C = h v iv Dh v + ... . (22) 

The collinear quark Lagrangian can be expanded as 

£ c = 4°) + £« + £f + .... (23) 

The first three terms are 

4 0) = L, P >{in-D+gn-A n>q +(f T +g4l q ) W±W* (fr+gflj)}^, 

4 2) = L{iprWlwHpr-(f T +g4l q )wlwl(n.iD)W^ (24) 

with = - igA£ s , where gives the order 0(A°) interactions [23 EH], and the expressions 
of C)P and of were derived in jHTJ and in \b'2\ . respectively. For the mixed effects, which are 
power-suppressed, the usoft-collinear Lagrangian has been derived up to 0(A 2 ) jHS]- Note that the 
results in [S3] represent an expansion of SCET in the hybrid momentum-position space. A manifestly 
gauge-invariant expansion in the position space has been derived in which each operator has a 
homogeneous power counting in A. 

After defining the power counting rules, I explain how to construct collinear factorization theorem at 
the operator and Lagrange level. The idea is to start with an operator relevant for a high-energy QCD 
process, which is characterized by some power of A. Draw the diagrams based on this operator and the 
effective Lagrangians in Eqs. (}2"2"|) and (|23p. The power of a diagram is the sum of those for the loop 
measures, propagators, vertices, and external lines. Note that powers in A shift between the propagators 
and vertices in different gauges. I will adopt the Feynman gauge the same as in Sec. 12.11 in order to 
compare the formalism in SCET and in perturbation theory. Those diagrams, whose contributions scale 
like the power the same as of the specified operator, contribute to the corresponding matrix element. 

Consider the diagrams in Fig. 0] The one-loop diagram in Fig. EJ^a) involves the leading-power 
operator for deeply inelastic scattering (DIS) [SB] . 

Odis = £ tn, P >W | C{P+, V-, Q, fi) Cn,p , (25) 
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\ 



(a) 



(b) 



Figure 4: Diagrams for comparing the power counting in SCET and in perturbation theory for the 
processes (a) DIS and (b) B — ► Dn, respectively. Dashed lines are collinear quarks, double solid lines 
are usoft or soft heavy quarks, and the single solid lines are soft light quarks. Gluons with a line through 
them are collinear, while those without a line are soft or usoft. 



with V± = ± V, and Q representing the momentum transfer from the virtual photon. The Wilson 
coefficient C, equivalent to the hard kernel, absorbs short-distance dynamics. When taking the proton 
matrix element of Odis 5 the V+ dependence of the Wilson coefficient leads to a convolution with parton 
distribution functions, which is the collinear factorization formula for DIS. Because of £ ~ A and 
W ~ A , the operator in Eq. (}25|) scales as A 2 . 

The two collinear gluon vertices come from the leading-power interactions of C^' in Eq. (|24|). Hence, 
Fig. Ufa) is characterized by the power, 



The factor outside the square brackets is for the external fields, the first term in the square brackets 
counts the collinear loop measure, and the second factor counts the three collinear propagators following 
Eq. (JBJ) or Table H The last factor in the bracket is the power of momentum in the quark-quark-gluon 
vertices in which are either (T, T) ~ (A, A) or (+, — ) ~ (A , A 2 ) in Feynman gauge (see the all-order 
proof of the collinear factorization theorem in Sec. 12. lj) . Equation (J26j) indicates that Fig. Ufa) is of 
the same power as the operator Odis, an d that nonperturbative collinear gluon exchanges of this type 
contribute to the leading-twist parton distribution function defined by (£n,p'W( , fi/2)W^ ni p) . 

It is easy to see that the above power counting is similar to that for Fig. Efb) in Sec. 12. II if one drops 
the power associated with the external collinear quark fields, Eq. (}2l)j) . being of 0(A°), corresponds to 
a logarithmic divergence, which should be absorbed into the distribution function. The strategies of 
the two approaches are compared below. In perturbation theory one starts with Feynman diagrams in 
full QCD. Look for the leading region of the loop momentum defined by Eq. in which one makes 
the power counting of the Feynman diagrams. It can be found that the approximate loop integral in 
the leading region is represented by a diagram of the type of Fig. Ufa) at 0(a s ). The distribution 
amplitude in Eq. (|T2*|) then collects this type of diagrams to all orders. In SCET one first constructs 
the various effective degrees of freedom describing infrared dynamics and the interactions in Eqs. (|22j) 
and (|23|). and defines their powers. Draw the diagrams based on the effective theory and then make 
the power counting. It can be shown that the diagrams of the type of Fig. Ufa) build up the leading- 
twist distribution amplitude in Eq. (|T2"j) . That is, one arrives at Fig. Ufa) through approximating loop 
integrals in the full theory in the former, but does at the operator and Lagrange level in the latter. 
Therefore, the derivations of collinear factorization theorem from both approaches are equivalent. 



At leading power, the external operator for the nonleptonic decay B — > Dtt in SCET is given by 




(26) 




(h ( :)ST h {1,T A 



}S j h^) $&WC {0JSi (P, P ] ) T e {1, T A } ^ , 



(27) 
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Figure 5: Collinear factorization of the B — > Dtx decay in SCET. 



where are the spin structures from the Fierz identity in Eq. ©. From Tabled 0{o,8} ~ A 5 is the 
base A-dimension for this process. In the three-loop diagram in Fig. Efb) all interactions are taken from 
the lowest-power Lagrangian 

4 0) - The direct power counting gives A : 

Here the first term counts the dimension of the external heavy quark fields and collinear quark fields. 
The term in curly brackets counts powers of A from the light soft spectator quark lines and the soft 
gluon propagator that does not participate in a loop. The factors in the first square bracket are the 
measure, propagators, and vertices for the soft loop. In the final square bracket the A factors are given 
for the measures, propagators, and vertices in the two collinear loops. The above power counting implies 
that Fig. IHb) contributes to the leading collinear factorization formula of the B — ► Dn decay. 

Note that collinear gluons do not attach the heavy quarks, which can not remain on-shell after 
emitting or absorbing collinear gluons [27|. Equivalently, no collinear divergence is associated with a 
massive particle in perturbation theory (SO]- Collinear gluons are not emitted by the soft spectator 
quarks in the B and D mesons either, since they do not produce pinched singularities (50! ■ The 
collinear divergences associated with the pion have been collected by the Wilson line W in Eq. (|27|) . 
Nonfactorizable soft gluons decouple from the pion due to the argument of color transparency [BT] 
mentioned in Sec 12.11 They contribute only at the subleading power. 

Following the above explanation, the proof of the collinear factorization theorem for the decay 
B — > Dn in SCET is trivial [66] : the leading-power diagrams involve soft gluons exchanged among the 
quarks in the B and D mesons, and collinear gluons exchanged between the quarks in the pion. The 
former give the B — > D form factor F BD , defined as the matrix element of the frist piece of 0{o,8}- 
The latter lead to the pion distribution amplitude <f> n (x), defined as the matrix element of the second 
piece with the Wilson coefficient C{ ,8} being excluded. The above discussion indicates that the pion 
distribution amplitude in Eq. (|18|) can be constructed in the framework of SCET. One simply identifies 
the correspondence of the quark fields q in Eq. (fTHj) and the collinear effective fields in Eq. (|27j). 
which are equivalent in the collinear region, and choose as 75 /n. The two collinear Wilson lines 
W correspond to the two pieces of path-ordered exponential in Eq. (|18|). The contributions from the 
diagrams in Fig. El are of the same power in A as of the effective current £jfywr t W*€Q in Eq. (|27|l . 

Hence, only the diagrams of the type shown in Fig. exist at leading power in SCET. The corre- 
sponding collinear factorization formula is then written as 

(D(*h\HjB) = NF B ^ DW (0) / dxT(x,E,fjt) <f>Jx,fi) + --- , (29) 

Jo 

where H w is the weak effective Hamiltonian, T(x, E, fi) a calculable hard kernel, and the ellipses denote 
terms that vanish faster than the leading term as the pion energy E — > 00. The dependences on 
the renormalization scale fi cancel between T and 4> n . The convolution in x is a consequence of the 
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determined by matching the effective theory to the lull theory. Equation 1)29)1 was proposed m [68J, 
proven to two- loops in jHH], and proven to all orders in a s in 66J. 

2.3 Jet Function 

The application of collinear factorization theorem to exclusive B meson decays has encountered a 
difficulty: the evaluation of the B — > 7r transition form factors suffers the singularities from the end 
point of a momentum fraction x — > [701 17JJ IHj • These singularities are logarithmic and linear in 
the leading-twist (twist-2) and next-to-leading-twist (twist-3) contributions, respectively. On the other 
hand, the double logarithms a s In 2 x from radiative corrections were observed in the semileptonic decay 
B — > nil/ [73] and in the radiative decay B — ► jlu [7j]. It has been argued that when the end- 
point region is important, these double logarithms are not small expansion parameters, and should be 
organized into a quark jet function systematically in order to improve perturbative expansion. The 
procedure is referred to as threshold resummation |75j . The resultant Sudakov factor is found to vanish 
quickly as x — > 0. It turns out that in a self-consistent perturbative evaluation of the B — > tt form 
factors, where the original factorization formula is convoluted with the jet function, the end-point 
singularities do not exist [To] . 

I take the radiative decay B — > --flu as an example. The momentum Pi of the B meson and the 
momentum P2 of the outgoing on-shell photon are parametrized as 

P 1 = ^|(l,l,O r ), P 2 =^(0,77,O T ), (30) 

where 7] denotes the energy fraction carried by the photon. Assume that the light spectator quark 
in the B meson carries the momentum k. Consider the kinematic region with small q 2 , q = P\ — P2 
being the lepton pair momentum, i.e., with large 7], where perturbative expansion is reliable. The four 
components of the spectator quark momentum k are of the same order as A. Here A represents a 
hadronic scale, such as the B meson and b quark mass diffetrence, mg- m&. In collinear factorization, 
only the plus comonent k + is relevant through the inner product k ■ P2. 

The lowest-order diagrams for the B — > 7/z/ decay are displayed in Fig. [TJ but with the upper quark 
(virtual photon) replaced by a b quark (W boson). It is easy to observe that Figs-H^a) and[T^b) scale like 
l/(Amfi) and l/m 2 B , respectively, implying that Fig.^b) is power-suppressed. Below I will concentrate 
on Fig. ^a). According to the leading-twist collinear factorization theorem discussed in Sec. 12.11 the 
B — > '-flu decay amplitude is written as the convolution of a hard kernel H(x,7],rriB, fi) with the B 
meson distribution amplitude + (x,/i) over the parton momentum fraction x = k + jP\ |76j . 

A(,, mB ) = Jd x M*,riH(*,,, mB ,,), (31) 

This expression has been derived in the framework of SCET [77J [7H], in which the hard kernel was 
further factorized into H = HhMi with the function Hh and Hi being characterized by the scale and 
i/Amj, respectively. Note that the jet function in [771 [THj, referred to ifj, differs from that in [75] . 

Equation ()3~T]) is appropriate for the region with k + ~ 0(A), in which the only infrared divergences 
are the soft ones absorbed into the B meson distribution amplitude [5UJ. Near the end point k + ~ 
0(A 2 /rriB), the internal quark in Fig. ^a) carries a large momentum P 2 — k with its invariant mass 
vanishing like (P 2 — k) 2 = —2xP\ ■ P 2 ~ 0(A 2 ). This kinematics is similar to the threshold region 
of DIS with the Bjorken variable xb — > 1, where the scattered quark also carries a large momentum 
and possesses a small invariant mass (1 — xb)s, s being the center-of-mass energy. In this region the 
scattered quark produces a jet of particles, to which the radiative corrections contain additional collinear 
divergences. Hence, a jet function needs to be introduced into the collinear factorization formula for DIS 
[79j . Similarly, a jet function has been incorporated into the factorization of direct photon production 
at a large photon transverse momentum (threshold) jHUJ- Here a jet function is associated with the 
internal quark near the end point of the momentum fraction involved in the decay B — > '-flu. 
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order correction to the weak decay vertex shown m rig. |2jdj. lhis divergence can be extracted by 
replacing the b quark line by an eikonal line in the direction of n. The factorization of the fermion flow 
is achieved by inserting the Fierz identity in Eq. (j3J), in which the first and last terms contribute in the 
combined structure Ujlik — > Iik(A A)ij/^- Assigning the identity matrix / to the trace for the hard 
kernel, one obtains Fig. ^a). The matrix jh jh/ 4 then leads to the loop integral [75] . 



\P 2 -k + If 



_nP_ 
n-lP 



= -^C F ln 2 x + --- (32) 

47T 

which are the same as those derived in [73J [THl [7S]. The correction to the photon vertex in Fig. [2{e) 
contains only the single logarithm a s In x, since the phase space of the loop momentum is restricted to 
< / + < k + ~ 0(A 2 /ms). The self-energy correction to the virtual light quark deos also. For the 
explicit expressions for the 0(a s ) corrections from Fig. [21 refer to |7fij . 

The all-order factorization of the jet function from the decay B — > 7/z/ has been proved following the 
procedure in Sec. 12. 1[ which provides a solid theoretical ground for the modified formalism appropriate 
for the end-point region. The jet function J(x) is defined via 



J(x)q(P 2 -k) = (q{P 2 - fc)|g(0)- jh fiexp 







ig / dzn ■ A(zn) 



|0) . (33) 



The spinor q{P 2 — k) is associated with the internal quark, through which the momeutm P 2 — k flows. 
It is then understood from Eq. ()33)1 that the jet function is universal. 

I then discuss threshold resummation of the double logarithms a s In 2 x in the covariant gauge d- A = 
0, which have been collected into the jet function to all orders. Threshold resummation for inclusive 
QCD processes has been studied intensively [EH IH3- Here I will adopt the framework developed in 
[HBl Ell, which has been shown to lead to the same results as in [HU [H2]- First, allow the vector n to 
contain a (small) minus component n~. This modification, regularizing the collinear pole, extracts the 
double logarithm as shown in Eq. (|32|). The definition in Eq. (j3*3*j) involves three variable vectors: the 
Wilson line direction n, the large momentum P 2 , and the spectator momentum k. The scale invariance 
in n, as indicated by the Feynman rule associated with the eikonal line along n, implies that the jet 
function must depend on k through the ratio n ■ k/n ■ P 2 . 

The next step is to derive the evolution of the jet function in x, i.e., in k + = xP^~ by considering 
the derivative, 

- dJ n - k .p«lL i (34) 



dk+ P 2 -k z dn a 

where the chain rule has been applied to relate the derivatives with respect to k and to n. The 
differentiation d/dn a operates on the eikonal line along n, giving 

n-k Da d n M A n-k P 2 ■ I 

P 2 -k 2 dn«n-l~ n-r % ~ P 2 ■ k n ■ I ^ ' 1 j 

The loop momentum I flowing through the special vertex does not generate a collinear divergence 
due to vanishing of the numerator P 2 ■ I in the special vertex n M . It is easy to confirm that the ultraviolet 
region of I does not produce lnx either. Therefore, one concentrates on the factorization of the soft 
gluon emitted from the special vertex, which can be achieved by applying the eikonal approximation to 
internal quark propagators, leading to n u /n- 1. Following the reasoning in [H3], the derivative of the jet 
function is written as 

dJjx) . 2n f dH g» u n v . + 

x — - — = —10 Up I -, — 7T — :J(x — I P\ ) , 36 

dx y J (2tt) 4 n-l I 2 n-l 1 1 ; ' v ; 

where the argument of J in the integral arises from the invariant mass of the internal quark, (P 2 —k+l) 2 w 
— 2(x — l + / Pi)P\ ■ P 2 . The integrand corresponds to the diagram with the soft gluon attaching the 
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equation, 



x- 



dJ(x) 



dx 



2n' 



c f 



d{ 



-AO , 



(37) 



where the variable change from l + to £ has been made. The plus distribution is defined such that, when 
l/(£ — x) + is integrated with a function /(£), one must replace it by /(£) — f(x) in the integral. It has 
been shown that the above evolution equation is similar to that for unintegrated parton distribution 
functions involved in inclusive QCD processes [85J, which resums the same double logarithm a s In 2 x. 
The analytical solution is a Sudakov factor, 



J{x) 



7T 

exp I — 7^ 




x 



,exp(t) 



sm 



7T 



-7 K t J exp ^- ^7«-*' 



(3* 



with the anomalous dimension 7^ = a s Cpjn. It is trivial to check that J(x) is normalized to unity, 
/ J(x)dx = J(l) = 1 86J. Obviously, Eq. (|3*8*|) vanishes at x — > 0, because the integrand is an odd 
function in t, and at x — > 1 due to the factor (1 — x) exp W. Moreover, Eq. ()38j) provides suppression 
near the end point x — > 0, which is stronger than any power of x. This is understood from vanishing of 
all the derivatives of Eq. (|38|) with respect to x at x — > jTS] . To the accuracy of the next-to-leading 
logarithms, the running of the coupling constant ot s should be taken into account, and Eq. (j3*S]l will be 
modified. However, the above features remain. 

I emphasize the differences among the Sudakov resummations for the B — > 7ZZ/ decay in the literature. 
In [23] it is the double logarithm In (hr/ms) that was resummed. In [25 it is the double logarithm 
In 2 (Ej/ tub), E 1 being the photon energy, that was resummed. In [ZH] the evolution from the scale of 
0(y/hmi) to the scale of 0(m&) was derived by solving the renormalizat ion-group equations, 



d 



dlnfi 



din fx 



-r cusp (a s ) ln^- 
r cusp (a s ) ln ^J 



Hh(fJ-) , 



(39) 



7 l + 



/•oo 

Hi(l+,fi)+ (kjT(u,l + ,a 8 )Hi(u,fi). (40) 
Jo 



r cusp is the universal cusp anomalous dimension familiar from the theory of the renormalizat ion of 
Wilson loops jHZj- The anomalous dimensions 7 and 7' are given by [ZH1 



7 (a.) = -2C F ^ + 0(a 2 s ) , 7 '(« s ) = 0(a 2 s ) . 
Arc 

The function T obeys J da; r(u;, a;', a s ) = 0, whose one-loop expression is written as 

'lu' e(uj-Lu') 9(u/-u)' 



(41) 



U) OJ — UJ 



UJ — UJ 



(42) 



It was found that the resummation effect decreases the magnitude of the radiative corrections, i.e., the 
renormalizat ion-group improved kernel is closer to the tree-level value than the one-loop result [75] . 

The formalism for threshold resummation has been extended to the semileptonic decay B — > ttIv in 
the fast recoil region of the pion. The B meson momentum Pi is the same as in the decay B — > '-flu, 
and the pion momentum P2 is the same as the photon momentum. Leading-twist factorization theorem 
for the B — ► 71 form factor F(q 2 ) has been proved in [oP] . 



F ( ( l 2 ) = Y\ I dxidx 2 (f> m (x 1 )H m {x 1 ,X2,r])(f) n (x2) 
Jo 



(43) 



which holds in the region with X\ ~ 0(A/m^) and with x<i 
amplitudes m will be defined in the next section. 



0(1). The light-cone B meson distribution 



14 



Figure 6: Leading-order contribution to F 



Since Fig. EJa), proportional to l/(xix%), is more singular at small x%, one considers the end-point 
region with X2 ~ 0(A/rriB), where the internal b quark propagator scales like l/(Ams). The loop 
correction to the weak vertex, where the radiative gluon attaches the virtual b quark and upper valence 
quark in the pion, generates the double logarithm a s In 2 X2 from the collinear region with the loop 
momentum parallel to This double logarithm, similar to that in Eq. (|32|) . is grouped into a jet 
funciton. It is easy to show that this jet function obeys the evolution equation in Eq. (J37j) . Hence, the 
threshold resummation leads to a result the same as Eq. (|38|). That is, the Sudakov factor is universal. 
The analysis for Fig. Efb) is similar to that for the decay B — > jlis. In the end-point region with 
x\ ~ 0(A 2 /m B ), additional collinear divergences associated with the internal light quark are produced. 
The loop correction to the weak vertex, where the radiative gluon attaches the b quark and the virtual 
light quark, gives the double logarithm a s ln 2 xi, whose factorization is the same as of Fig. 12(d). 

The modified collinear factorization formula appropriate for the end-point region is then written as 

F (q 2 ) = EE/ dx 1 dx2(f)m{xi)H^(xi,X2,r])J{xi)4 > n{x2) , (44) 

m=+- i=l,2 J ° 

with the index % — 1 (2) corresponding to Fig. Efb) [(a)]. If J(x) is excluded, the above expression is 
divergent because of H$ (x 1/x 2 . and oc x 2 at x 2 — > 0. Including the threshold resummation, the 
B — > n form factor is calculable without introducing any infrared cutoffs [7U1IZ2]- The numerical effect 
from the jet function on the B — > n form factor has been examined in [To] . 

In a recent work based on SCET, a jet function has also been defined in the analysis of the decay 
B — > nlv jSHl- It was also concluded that the end-point singularity does not exist in the B — > n transiton 
form factors in the convolution with the jet function. I stress that the jet function in .88j differs from 
the one considered here, and that the smearing mechanism of the end-point singularity is also different: 
it is not attributed to the Sudakov mechanism discussed above. The jet function in jSH], absorbing 
dynamics characterized by 0(y/ Am B ), more or less corresponds to the finite piece of the hard kernels 
in collinear factorization theorem without threshold resummation. In the case of the B — > n transiton 
form factors, it can be identified as the piece from Fig. EJa), which is proportional to l/a^- This piece 
is free of the end-point singularity, nothing to do with the Sudakov effect. This point will be elucidated 
in detail in Sec. 14.31 

3 kr Factorization 

Both collinear and kj- factorizations are the fundamental tools of QCD perturbation theory, where 
kr denotes parton transverse momenta. For inclusive processes, consider DIS of a hadron, carrying 
a momentum p, by a virtual photon, carrying a momentum q. Collinear factorization [89J and kx 
factorization |§U1 WI\ apply, when DIS is measured at a large and small Bjorken variable xb = 
—q 2 /(2p ■ q), respectively. The cross section is written as the convolution of a hard subprocess with a 
hadron distribution function in a parton momentum fraction x in the former, and in both x and kr 
in the latter. When x B is small, x > xb can reach a small value, at which k? is of the same order 
of magnitude as the longitudinal momentum xp, and not negligible. For exclusive processes, such as 
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previous section, lhe range ol a parton momentum traction x, contrary to that m the inclusive case, 
is not experimentally controllable, and must be integrated over between and 1. Hence, the end-point 
region with a small x is not avoidable. If no end-point singularity is developed, collinear factorization 
works. If such a singularity occurs, indicating the breakdown of collinear factorization, kp factorization 
should be employed. In fact, the recent observation QF 2 (Q 2 ) / F^Q 2 ) ~ const. -Pi and F 2 being 
the proton Dirac and Pauli form factors, respectively, indicates that k? factorization is the appropriate 
tool for studying exclusive processes jM] ■ Since k T factorization theorem was proposed jHSl EH] , there 
had been wide applications to various processes |§T] . 

In this section I review kp factorization theorem for exclusive processes. It is more convenient to 
perform kp factorization in the impact parameter b space, in which infrared divergences in radiative 
corrections can be extracted from parton-level diagrams explicitly. The procedure is basically similar 
to that for collinear factorization in Sec. 12.11 if the proof is performed in the impact parameter b space. 
It has been observed that collinear factorization is the b — > limit of k? factorization. I explain how to 
construct a gauge-invariant 6-dependent meson wave function defined as a nonlocal matrix element with 
a special path for the Wilson line. The application of hp factorization theorem to exclusive B meson 
decays, and the behavior of 6-dependent B meson wave functions are discussed. Retaining the parton 
transverse degrees of freedom, the double logarithms a s In 2 kr appear, which should be organized to 
all orders. The basic idea for kp resummation of these double logarithms into a Sudakov factor [21] 
is given. The end-point singularity in the heavy-to-light transition form factors can also be smeared 
[21 123 H3 EHl El] by including this Sudakov factor. 



3.1 Gauge Invariance 

I again start with the process rcy* — > 7 [100 . This process, though containing no end-point singularity, is 
simple and appropriate for a demonstration. The momentum Pi (P2) of the initial-state pion (final-state 
photon) is chosen as in Eq. (fl|). I explain how to perform the factorization of the collinear enhancement 
from I parallel to P\ without integrating out the transverse components It- The lowest-order diagrams 
are displayed in Fig.Q and the 0(a°) kr factorization formula is the same as the collinear factorization 
formula in Eq. P|). That is, none of ft°\ and depends on a transverse momentum. The wave 
function and the hard kernel become /r-dependent through collinear gluon exchanges at higher orders. 

The 0(a s ) kp factorization formula is a sum over the diagrams in Fig. |21 the same as Eq. (J7J), 
but with each term being written as the convolution in the momentum fraction £ and in the impact 
parameter b, 

= j ' d^-^-^f\x^b)H^{^b) + ft°\x)H?\x) . (45) 

The above expression, with the 0{a s ) wave functions §f\x, £, b) and H^(^,b) specified, defines the 
0(a s ) hard kernels Ti^ix), which do not contain collinear divergences. Equation (j4T)|) is a consequence of 
the assertion that partons acquire transverse degrees of freedom through collinear gluon exchanges: H.^ , 
convoluted with the lowest-order /^-independent ft ', is then identical to that in collinear factorization. 
As shown later, this consequence is crucial for constructing gauge-invariant hard kernels. 

The 0(a s ) wave functions obtained from Figs.|2fa) and|2fc) are the same as in collinear factorization. 
The kr factorization of Fig. [2](b) leads to the wave function, 




(46) 

The Fourier transformation introduces the additional factor exp(— i\? ■ b) into the wave function $j 
compared to the result in collinear factorization in Eq. (jHJ), since the hard kernel depends on in this 
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Figure 7: The special path for the Wilson line in a 6-dependent wave function. 



case. The 0(a s ) wave function extracted from Fig. 12(d) is written as 

ig 2 C F r dH 5 x pi+ / 1 n" 



4P X + 7 (2tt) 4 ^ y ' ' (xPi + 2 'l 2 n-l 

(47) 



- *) - 5 - x + lp J e -^- b 



The second term acquires the additional factor exp(— Ht ■ b) from the Fourier transformation, because 
it corresponds to the case with the loop momentum I flowing through the hard kernel. It is easy to 
observe that the soft divergences cancel among the 0(a s ) radiative corrections: in the soft region of I we 
have exp(— ily • b) pa 1 and l + pa 0, and the two terms in Eq. ()47j) cancel. Similarly, the soft divergences 
cancel among Figs. 2(a)-2(c). 

One constructs the parton-level wave function as the nonlocal matrix element in the b space, 



$(x,£,b) = i I ^-e-^v" (0\q(y) l5 fiPexp \-ig f ds ■ A(t 
J 2ir L Jo 



q{0)\q{xP 1 )q(xP 1 )) , (48) 



with the coordinate y = (0, y~,h). The path for the Wilson line is composed of three pieces: from to 
oo along the direction of n, from oo to oo + b, and from oo + b back to y along the direction of — n 
as displayed in Fig. [7| The first (third) piece corresponds to the eikonal line associated with the first 
(second) term in Eq. (}4Tj) . 

For the evaluation of the lowest-order hard kernel, one neglects only the minus component l~ in the 
denominator [see the second term on the right-hand side of Eq. (jl(Jj) ]. 

(P 2 _ xPl + if « _( 2 ^ • p 2 + ll) . (49) 

Note that in collinear factorization both l~ and It are dropped. The 6-dependent hard amplitude is 
then given by, 



#"»(?,'>) = />iT« (0, (?,ir)exp(a r .b), 



mtM = ^ h ^p f? 5> • (50) 

Equivalently, the above H.( \£,,It) is derived by considering an off-shell q quark, which carries the 
momentum £P X — It, and the leading structure Pi 75 associated with the pion, which is the same as in 
collinear factorization. 

I now demonstrate the gauge invariance of fc^ factorization theorem. Equation (}4*Sj) is explicitly 
gauge-invariant because of the presence of the Wilson link from to y [SHHH]. H {1) (x), the same as 
in collinear factorization, is gauge-invariant. From Eq. ()45jl . the gauge invariance of $ (1) (x,£,6) stated 
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H ( >{!;, b). similarly, the kj> iactonzation iormula of U{a s ), 

g^\x) = J d^^y 2 [^ 2 \x,^bM \^b) + ^ 1 \x,^bM 1 \^b)]+^\x)n (2 \x) , (51) 

leads to the gauge invariance of H^(£,b). Therefore, the hard kernels in hp factorization are gauge- 
invariant at all orders. 

After determining the gauge-invariant infrared-finite hard kernel H, one convolutes it with the 
physical two-parton pion wave function, whose all-order gauge-invariant definition is written as 

%(x,b,Q,pL)=i [ ^- e - ixP ^~ (0\q(y) l5 fiV exp \-ig f ds ■ A(s)} q(0)\*(Pi)) ■ (52) 

J Z7T L JO 

The relevant form factor F(Q 2 ) for the process 717* — > 7 is then expressed as 

F(Q 2 ) =£dxj 7^4$ir(a:, b, Q, a)H(x, 6, Q, u) . (53) 

where both the dependences on Q and on the factorization scale \x have been made explicit. It has been 
concluded that predictions derived from kr factorization theorem are gauge-invariant and infrared-finite 

una. 

In summary, a two-parton 6-dependent wave function is factorized from parton-level diagrams in a 
way the same as in collinear factorization (for example, under the same eikonal approximation), but the 
loop integrand is associated with an additional Fourier factor exp(— Ht ■ b), when the loop momentum 
I flows through a hard kernel. A fcr-dependent hard kernel is obtained in a way the same as in collinear 
factorization, but considering off-shell external partons, which carry the fractional momenta k = xP—hx 
(k 2 = —k T ), P being the external meson momenta. Then Fourier transform this hard kernel into the b 
space. The insertion of the Fierz identity to separate the fermion flow between a wave function and a 
hard amplitude is the same as in collinear factorization. For inclusive processes in small xb physics, the 
gauge invariance of the unintegrated gluon distribution function and of the hard subprocess of reggeized 
gluons, being also off-shell by —k 2 -,, is ensured in a similar way. The distinction is that the structures of 
7-matrices from the Fierz identity are replaced by eikonal vertices, which contain only the longitudinal 
components p?!] . 

3.2 B Meson Wave Functions 

In this subsection I review the k^ factorization theorem for exclusive B meson decays by considering the 
radiative decay B — ► '-flu. It has been shown that in heavy quark limit a gauge-invariant 6-dependent 
B meson wave function can be defined, which absorbs soft divergences in the decay process, differing 
from the collinear divergences in the pion wave function. As explained below, exclusive B meson decays 
are characterized by the scale i/Am^. In terms of the power counting in SCET, the soft dynamics 
discussed here is referred to as the usoft one, since the typical momentum behaves like [SU] 

P~(A,A,A)~m B (A 2 ,A 2 ,A 2 ), (54) 



for the expansion parameter A ~ yA/m^. It is possible to construct a light-cone B meson wave 
function, if an appropriate frame with the photon moving along the light cone is chosen. 
Figure Hl^a) gives the parton-level amplitude, 

gl°\x) = eq(k) / ( ^ 2 ~^ 2 7m(1~ 75)^1 -fe), (55) 

which does not depend on a transverse momentum. Inserting the Fierz identity in Eq. Q into Eq. ([55)1 . 
one obtains Eq. (0J with 

H (o ){x) = MA ^(i-i 5 )(A + m B )(A/V2h 5 ] ^ (56) 

2xP± ■ P2 
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leading- twist wave iunctions associated with the structures 757 . ror the B — > 7/1/ decay, only the 
structure 757" 1 " = 75 /n contributes: since /e in Eq. (jS^jl involves 7t, only the structure 7~75 = ^75 
contributes to the hard kernel. 

Next one considers the 0(a s ) radiative corrections to Fig. da) shown in Fig. |21 and discuss the 
factorization of the soft divergence from the region of the loop momentum in Eq. (|54j) . The dependence 
of the B meson wave function on the transverse momentum is generated by soft gluon exchanges. The 
analysis is similar to that in Sec. 13.11 and one derives Eq. (J45|) . The factorization of the two-particle 
reducible diagrams in Fig. |2t&)-G(c) is straightforward. Take Fig. EJb) as an example. Employing the 
eikonal approximation in the heavy quark limit, one has 

^^^"S^-^^-^ (57) 



(Pi — k + I) 2 — m 2 v-l 
with the velocity v = Pi/ tub- The 0(a s ) wave function extracted from Fig. 12(b) is then written as 



•*"(*.«.») = V/^^t^OV*^-*^^-^^ 6 ■ (58) 

Performing the contour integration over, say, l~, one observes that the integral is singular only when 
the component Z + is of 0(A). This observation implies that the infrared divergence associated with the 
B meson is of the soft type, and that / + , being of the same order of magnitude as k + = xPi~, is not 
negligible in the ^-function. Therefore, the soft divergences in Figs. E(a)-|2(c) do not cancel in B meson 
decays [24J. The explanation is simple: the light spectator quark, carrying a small amount of momenta, 
forms a color cloud around the b quark. This cloud is also huge in space-time, such that soft gluons 
resolve the color structure of the B meson. 

Diagrams with the radiative gluon attaching the internal quark in Figs. 12(d) and |2(e) also contain 
soft divergences, because the internal quark is off-shell by O(Am^), which defines the characteristic 
scale of the decay B — ► 7/z/. Note that the internal quark in the process 7r7* — > 7 is off-shell by 0(Q 2 ). 
One extracts the 0(a s ) wave function from Fig. 12(d), 

(!) —ig 2 CF f d A l _. . . _ , 1 n-v 



*SW.6) - ^/^ l)75W -^ 2n , w , 




S(Z-x)-5[Z-x + — )e-^ h 



(59) 



The eikonal approximation in Eq. (|57|) has been applied. The above expression implies that the infrared 
divergences in the irreducible diagrams can also be collected by the eikonal line along the light cone. 
This is attributed to the choice of the frame, in which the photon moves in the minus direction. 

Following the procedure in Sec. 13.11 one constructs a gauge-invariant light-cone B meson wave 
function, 



$ + (x,6,m B ,/i) = i [^ e - ixP i y -{0\q(y)w + Vexv\-ig f ds ■ A{s) 

J Z7T JO 



b v (0)\B(Pi)) , (60) 



where b v is the rescaled b quark field, and the decay constant fs has been omitted. The Feynman rules 
associated with b v are those for an eikonal line in the direction of v defined in Eq. ([57]). The lowest-order 
hard kernel in the b space is given by Eq. (|50j) with 

wow , \ Mt /P 2 7M(l-7 5 )(/Pi + m g )(AV2)7 5 ] 

H (Ut) = ~ 6 2^ ■ P 2 + l\ ' (61) 

^ = (k + — l + ) / Pi being the momentum fraction. The above expression can be derived by considering 
an off-shell q quark of the momentum (£Pi~, 0, —It), and by contracting the parton-level diagram with 
the leading structure (/Pi + m^)(/i/v / 2)7 5 , which is the same as in collinear factorization. 
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whose all-order prooi can be performed m a similar way. ror this mode, both the leadmg-twsit B meson 
wave functions <j)±, associated with the structures 757 ± , contribute [50J. Moreover, contributions from 
the pseudo-scalar and pseudo-tensor two-parton twist-3 pion distribution amplitudes are also leading- 
power [HE]- The factorization of the corresponding collinear divergences has been proved [301 - The point 
is to replace the Dirac structure 757" by the corresponding ones 75 and 75 c™^ in Eq. (j3J). 

I then discuss the behavior of the B meson wave functions constructed in Eq. (J60j) . In the heavy 
quark limit the two-parton light-cone wave functions z 2 ) are defined in terms of the nonlocal 

matrix element [103, 104 : 



(0|?(y)r&„(0)|fl(Pi)> = -'^±tr 



2t 



(62) 



with t = v ■ y, y 2 = —b 2 , and T being a Dirac matrix . 

Consider the light-cone distribution amplitudes in terms of the variable uj = xnis [103J, 

</> ± (u) = lim ® ± (uj,y 2 ) , (63) 

where the wave functions $±(w, y 2 ), defined in Eq. (|60j) . come from the Fourier transformation of 
&±(t,y 2 )- The differential equations are written as |l()5j 

(u - 2A) <f> + (uj) + ujcp^uj) = J(uj) , (64) 
where and J(uj) denote the source terms due to three-parton wave functions ^/a, an d Xa' 

A fUJ roo AC f) 

The solution can be decomposed into two pieces: 

± H = 0f ) H + 0?M. (66) 

The functions (f>±(ui) are the solution with I(u>) = J(oj) = 0, corresponding to the "Wandzura-Wilczek 
approximation" |1()6^ ll()7j = = Xa = 0. The functions 4>±\u) are induced by the source terms 
I(uj) and J(u). The analytic expressions for the Wandzura-Wilczek part are given by 

4"V)- A± ^r A W (2A-o,). (67) 

The expressions for <f>±\ in terms of \E% and Xa, can be found in |108j . Equation (jBTjl is quite 
different from the model distribution amplitudes appearing in the literature. One example of such 
models is 11031 





exp 


f UJ 


UJq 




\ ^0 






< UJ 


4> G - N (") = - 


exp | 




uj 




v ^0 



(68) 

with uj q = 2A/3, which are inspired by the QCD sum rule estimates 103J. Note that, however, the 
behavior <p9 N (u) ~ uj and (f)_(uj) ~ constant at uj — ► is consistent with Eq. ([67)1 . Another example 
comes from solving an integro-differential equation |109j : evolution effects generate a radiative tail, 
which falls off slower than 1/uj. 
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f AN 9 3 

(«-A)£ + 2 



-$_)=0, 

<9cj 2 <9;Z 2 

d 2 d$ + 
_$_ + 2 

2 ctaj 2 cfe 2 



$_) +2 



d 2 



So; 2 \ cfe 2 



dz 2 



(69) 



The solution in the Wandzura-Wilczek approximation is given by 



w,6) 



>/6 



cos 



^/u;(2A 



6- 



7T 



(70) 



where the 6 dependence corresponds to 5{k^ — u;(2A — u)) in /c T space. It is observed that the longi- 
tudinal and transverse momentum dependences do not seperate (factorize) in Eq. (|7D|L contrary to the 
assumption in many models [25 | I104[ TllOj. The Gaussian distribution for the fcr-dependence has been 
adopted in [2*5] . which exhibits strong damping at large b as exp (~u 2 3 b 2 /2) (see Sec. 14. Hj) . In contrast, 
the results in Eq. (|7()jl show slow-damping with oscillatory behavior. I mention that, despite of the 
different functional forms, the numerical results of the B — > it form factor derived from the B meson 
wave function in [53] and from Eq. (|70|) are very similar 



3.3 kx Resummation 

The inclusion of parton transverse degrees of freedom introduces a soft logarithm a s In b. Its overlap 
with the original collinear logarithm leads to a double logarithm a s \n 2 (Qb). This large logarithm must 
be organized in order not to spoil perturbative expansion. I explain the idea of kx resummation by 
taking the pion wave function as an example. It is known that single logarithms can be summed to all 
orders using renormalization group methods, while double logarithms are organized by the technique 
developed in |112| IllHj . I choose the axial gauge n ■ A = 0, in which the two-particle reducible dia- 
grams, like Figs.|2(a)-|2fc), contain the double logarithms, while the two-particle irreducible corrections, 
like Figs. Efd) andEfe), contain only single soft logarithms. If the double logarithms appear in an 
exponential form & w ~ exp [—const, x \n 2 (Qb)], the task will be simplified by studying the derivative 
of <& n , d^^/dlnQ = C§ n . It is obvious that the coefficient C contains only large single logarithms, 
and can be treated by renormalization group methods. Therefore, working with C one reduces the 
double-logarithm problem to a single-logarithm one. 

Consider the pion wave function <& n (x, b, Q, fx) defined in Eq. (J52j) . The two invariants appearing in 
Q n are Pi ■ n and n 2 , where n is allowed to vary away from the light cone. By the scale invariance of n 
in the gluon propagator, 



n • I (n ■ l)' 



(71) 



§ n depends only on a single large scale v 2 = (Pi ■ n) 2 /n 2 . It is then easy to show that the differential 
operator d/dlnQ can be replaced by d/dn: 

d n 2 „ d , s 

= ■ (72) 



dlnQ " Pi-n 1 dn c 

The motivation for this replacement is that the momentum Pi flows through both quark and gluon 
lines, but n appears only on gluon lines. The analysis then becomes simpler by studying the n, instead 
of Pi, dependence. 
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d -N" v = -J—(Ni"*l v + N va F) . (73) 



dn a n ■ I 

The momentum I that appears at both ends of a gluon line is contracted with the vertex, where the 
gluon attaches. After adding together all diagrams with different differentiated gluon lines and using 
the Ward identity in Eq. (|13|). one arrives at the differential equation of and the result contains 
the special vertex |2"4"j . 

9 P - j Pla • (74) 

An important feature of this special vertex is that the gluon momentum I does not lead to collinear 
enhancements because of the nonvanishing n 2 . The leading regions of I are then soft and ultraviolet, in 
which the subdiagram containing the special vertex can be factorized from the new function The 
left-over part is exactly Q^, and the subdiagram is assigned to the coefficient C introduced before. 

Therefore, one needs a function K to organize the soft divergences and G to organize the ultraviolet 
divergences in the subdiagrams. The differential equation of <$> n is then written as, 



b, Q, fj.) 



dlnQ 



2K(bfi) + - G{xv/pL) + - G((l - x)u/fi) 



(x, b, Q, fj) . (75) 



The functions K and G have been calculated to one loop, and the single logarithms have been organized 
to give their evolutions in b and Q, respectively jHS]- They possess individual ultraviolet poles, but their 
sum K + G/2 is finite, such that Sudakov logarithms are renormalization-group invariant. Substituting 
the expressions for K and G into Eq. (jZEJ), one derives the solution, 



® n (x, b, Q,n) = exp 



E 8(&,t>,Q) 



$«{x,b,n), (76) 



where the initial condition of the Sudakov evolution, <$> n (x, b, /x), contains the single- logarithm evolution 
in /i, and the intrinsic dependence on b |114j . The distribution amplitude (j) n (x, /x), defined in Eq. (fTHj) . 
is the 6^0 limit of $ n (x, b, //). The explicit expression for the exponent s, grouping the double 
logarithms, is referred to |25j . 

Note that the vector n has been varied away from the light cone in the above technique. The 
leading- logarithm resummation, being independent of the n, is still gauge invariant. The n dependence 
indeed appears in the next-to- leading- logarithm resummation for the wave function |95 | I104| . such that 
this piece becomes gauge dependent. However, this n dependence will be cancelled by that from the 
resummation of nonfactorizable soft gluons, which is also next-to- leading- logarithm |%3*| 1115] . That is, 
in a complete next-to- leading- logarithm resummation, the Sudakov factor is gauge invariant. 

Variation of exp(— s) with b and x is displayed in Fig. |H1 which shows a strong falloff in the large b 
and large x region, and vanishes for b > 1/Aqcd- Hence, Sudakov suppression selects components of the 
pion wave functions with small spatial extent b, and makes the hard scattering more perturbative. Once 
the main contributions to the factorization formula come from the small b, or short- distance, region, 
perturbation theory becomes relatively self-consistent. 

The above formalism has been generalized to the B meson wave function. In the axial gauge only 
the two-particle reducible diagrams generate the double logarithms. Figure a), giving the self-energy 
correction to the massive b quark, produces only soft enhancement, and is subleading. If the component 
k + of the spectator momentum is as small as 0(A), collinear divergences in Figs. Efb) andEfc), which 
arise from the loop momentum with a large component parallel to k, will not be pinched, and they also 
give only soft enhancements. This is consistent with the physical picture that the soft light quark can 
not interact with the heavy quark through a fast moving gluon. If there is nonvanishing probability of 
finding the light spectator with k + being of O(mg), such as in the model with a power- law decrease in 
k + , Figs. [21(b) and [21(c) contribute large double logarithms. Most of the models for the B meson 
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Figure 8: The Sudakov factor exp[— s(x, b, Q)] for Aq CD = 250 MeV. 



wave function in the literature favor k + ~ 0(A). That is, the k^ resummation for the B meson is 
not important. However, I will discuss this resummation, and allow the behavior of the B meson wave 
function to determine whether its effect is essential. 

The major difficulity in summing up the double logarithms in Figs. Efb) and |2fc) arises from the 
many invariants that can be constructed from Pi, k and n, such as P 2 , P± ■ k, Pi ■ n, k ■ n and n 2 . 
In the pion case the invariants are only Pi ■ n and n 2 . The fact that the B meson wave function 
contains many invariants fails the replacement of d/dk + by d/dn, because some large scales like P 2 
can not be related to n. Fortunately, this difficulity can be overcome by applying the heavy quark 
approximation in Eq. (JSTj) . This approximation also holds for collinear gluons with momenta parallel 
to k, since collinear divergences are independent of the direction of the eikonal line that collects the 
collinear gluons. Different directions correspond to different shifts of finite contributions between the 
wave function and hard kernels, i.e., to different factorization schemes |115j . However, it was argued 
|104j that the approximation in Eq. (J57|) is not suitable for collecting collinear gluons. 

Substituting the eikonal line along v for the b quark line, self-energy diagrams like Fig. |21^a) are 
excluded by definition |117j . The eikonal approximation also reduces the number of large invariants 
involved in the B meson wave function. We have the scale invariance in Pi in addition to the scale 
invariance in n. Hence, Pi does not lead to a large scale, and the only large scale is k + , which must 
appear through the ratios (k ■ n) 2 /n 2 and (k ■ v) 2 /v 2 . At leading-logarithm accuracy, the second scale 
does not appear. This observation can be verified by evaluating the soft function K and the hard 
function G for $ + However, the above argument for the survival of a single large scale has been 
questioned jlU4j . The Sudakov effect associated with the B meson is not important, and the dispute 
does not affect the numerics discussed in the following sections. 

Since $ + depends only on the single large scale (k ■ n) 2 /n 2 , the derivation reduces to the one in 
analogy with the pion case. One obtains 

$+(x, b, m B , p) = exp [s(x, b, m B )\ ®+(x, b, fi) , (77) 

with the same exponent s. The intrinsic b dependence, which is more important for a heavy meson, has 
been included into the initial condition of the Sudakov evolution, $+(x, b, fi). The behavior of $+(x, b), 
ignoring the single-logarithm evolution in /i, has been discussed in the previous subsection. 

4 Semileptonic and Radiative Decays 

The B meson decay constant and transition form factors, involving the hadronic effects in semileptonic 
and radiative decays, provide the nonperturbative inputs of many QCD methods. In this section I 
review the recent studies of these topics in LCSR, lattice QCD, PQCD QCDF, SCET, and LFQCD. I 
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quark mass can be determined by means oi a two-point correlation tunction similar to that lor the heavy 
meson decay constant. On the other hand, the above quantities are not very relevant to the leading- 
power formalism of nonleptonic B meson decays discussed in Sec. 5. The results for corresponding D 
meson decays will be quoted for comparison. 



4.1 Light- Cone Sum Rules 

QCD sum rules |118t I119j have been applied to various problems in heavy flavor physics. The idea is 
to calculate a quark-current correlation function and to relate it to hadronic parameters via dispersion 
relations. Take the B meson decay constant /# as an example |12U[ I121[ I122j , which is defined via the 
matrix element (0|m^gi7 5 6|i?) = /gWg, q = u,d. Consider the correlation function of two heavy-light 
currents, 

n(g 2 ) = z J d 4 ye iq - y (0\T[m B qi l5 b(y),m B bi l5 q(0)]\0) . (78) 

The amplitude II(g 2 ) can be treated by operator product expansion (OPE) at the quark level, if q 2 is 
far below m 2 , or parametrized as a sum over hadronic states including the ground-state B meson for 
q 2 > m 2 B . Assuming the quark-hadron duality, the expressions in the above two regions are related. 
Therefore, on the left-hand side of the sum rule, one has 

n(g 2 ) = -#^ + --- , (79) 

m B — q 2 

where the contribution of the ground-state B meson has been singled out, and • • • represents those from 
the excited resonances and from the continuum of hadronic states with the B meson quantum numbers. 
On the right-hand side of the sum rule, we have the expansion including the perturbative series in a s 
and the quark, gluon and quark-gluon condensates. A simple explanation of the quark-hadron duality 
has been given in |123t 1124] . Inserting the values of a s , and the condensates (G 2 ) and (qq) into the 
above sum rule, one estimates 

LCSR |125j . employed frequently for studying exclusive B meson decays, is a simplified version of 
QCD sum rules. Consider the B -> vr transition form factors [ESI Q2E1 H23 E2E1 Q22I HSDj, for which 
the correlation function is chosen as 



i jd A ye^ y (n + (P 2 )\T[u lx b(y),m b h l5 d(0)}\0) 

= F((P 2 + q) 2 , q 2 )P 2 , + F((P 2 + q) 2 , q 2 )q, . (80) 

Compared to Eq. (J78|) . the final state has been specified as a pion, and the twist expansion has been 
applied. The presence of the heavy quark mass justifies the twist expansion. 

At large virtuality | (P 2 + q) 2 — m 2 |3> Aq CD and q 2 <C m 2 , the correlation function is treated by OPE 
near the light-cone y 2 = 0. The perturbative part involves a convolution with the pion distribution 
amplitude 4> n (x) according to collinear factorization theorem in Sec. 12.11 The evaluation becomes 
simpler: it contains an integral only over the one-dimensional momentum fraction x, instead of over 
the four-dimensional loop momentum. The price to pay is that higher-twist contributions need to be 
included in terms of inverse powers of (P 2 + q) 2 — m 2 . On the hadron side, one has 

F((P I n\ 2 n 2 \- 2 f BF +^ m B i (o-i\ 

F((P2 + g) > q) ~ m 2 B -(P 2 + q) 2+ -- ' (81) 

where the ground-state contribution from the B meson contains a product of fs and the B — > tt form 
factor F + (q 2 ). The form factor F + , along with another one Fq, are defined via 



(n + {P 2 )\q 1 ,b\B{P l )) = F + {q 2 



2 2 

m n — ml 
(Pi + P2), ~ B q2 



+ F^H^UhL (82) 
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Ihe resulting sum rule is written as 



exp 




F (2 V,M 2 ,m fe 2 , So B ,/i) 

(83) 



2m 2 B f B 

±^if V, M\ ml si p) + FgWtf, M\ mg, sg, 

The mass scale M is associated with a Borel transformation usually performed in sum rule calculations. 
The scale /i is the factorization scale separating soft and hard dynamics. The effective threshold 
parameter Sq sets the lower bound of the B meson invariant, {Pi ± q) 2 > s B , above which the quark- 
hadron duality is assumed to hold. Long-distance dynamics characterized by scales lower than \x is 
absorbed into the universal nonperturbative pion distribution amplitudes. The first two terms on the 
right-hand side of Eq. ()83|) represent the twist-2 contributions up to next-to-leading order. The third 
term represents the leading-order twist-3 and twist-4 contributions. 

(2) 

For illustration, the leading term F is given by 



i 2 \q\ Ml ml s$, n) = m\U J ^ exp (- 

A \ 



Ff>( q lMlml,sln) = miU / -exp ( - m » ^ x) ) . (84) 



The lower integration boundary A = {m 2 — q 2 )/{sQ — q 2 ) originates from the subtraction of excited 
resonances and continuum states from both sides of the sum rule, which contribute to the dispersion 
integral in the B meson channel. These contributions are identical on both sides of the sum rule because 
of the quark-hadron duality assumed above. The explicit expressions for the remaining terms and 
Fo 3 ' 4 "* can be found in |126| I127j and in |128t I129j , respectively. The radiative correction to the twist-3 
cobtribution has been available |131j . 

The D-meson decay constant fo can be derived by a simple h —> c (B —>■ D) replacement in the 
sum rule for fg in Eq. ([79)1 . and by the necessary adjustment of the renormalization scale. One can also 
predict the ratios /b s //b an d /d s //d by setting the quark field q = s in Eq. (J7HJ). The values of fg and 
fo are sensitive to the b and c quark pole masses, m& and m c . Varying these masses in the intervals, 



one obtains |132j 



m b = 4.8 ± 0.1 GeV, m c = 1.3 ± 0.1 GeV, (85) 



/b = 170 =f 30 MeV, f D = 180 =f 30 MeV , 
f Bs /f B = 1-16 ±0.09, f D Jf D = 1.19 ±0.08. (86) 



Within uncertainties, the predictions in Eq. (J86|) agree with the lattice determinations of the heavy 
meson decay constants quoted in the next subsection. 

The LCSR predictions for F + |133j are presented in Fig. El This calculation includes twist-2 (leading- 
order and next-to-leading-order) and twist-3, 4 effects. The twist-2 and twist-3 contributions are roughly 
equal. The twist-4 contribution is less than 10% in the fast recoil region. The results are insensitive to 
the nonasymptotic behavior of the pion distribution amplitudes. At the maximal recoil, one finds [133 

Ff 7r (0) = 0.28 ±0.05 , Ff^O) =0.65 ±0.11 . (87) 

Note that QCD sum rules have a limited accuracy due to the truncation in OPE, to the duality 
approximation, to the variation of the corresponding auxiliary parameters, such as the Borel mass 
M, and to the contributions of excited states. A detailed discussion on the uncertainty from the above 
sources can be found in |133j . Moreover, large radiative correction to the B meson vertex, which reaches 
35% of the full contribution, or about half of the leading-order contribution, has been noticed in the 
correlation function in Eq. (J80j) . This 0(a s ) correction renders the sum rule for f B F + quite unstable 
relative to the variation of input parameters |126| I130j . This is the reason one considers the sum rule for 
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Figure 9: The LCSR prediction for the B — > n form factor F + (p 2 ) [133]. The full curves indicate the 
theoretical uncertainty, and the points represent various lattice QCD calculations from FNAL |134j (full 
circles), UKQCD [T33| (triangles), APE USHj (full square), JLQCD [H7j (open circles), and ELC 
(semi-full circle). 



fs at the same time in order to stabilize the sum rule for JbF+- the sum rule for also receives large 
radiative correction to the B meson vertex, such that the two large vertex corrections cancel in the ratio 
fsF+l fs |138j . However, the radiative correction to fs then becomes large. Therefore, an evaluation 
of 0(a 2 ) corrections to both the sum rules is necessary. Progress has been made in the calculation of 
the three-loop radiative corrections to the heavy-to-light correlator 139J. 

Replacing the pion with the kaon (including the distribution amplitudes and the decay constants) 
in the correlation function in Eq. (JHUj) . one obtains LCSR for the B —> K form factor, and the ratios 
[TUB] . 
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Bn 
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where the uncertainty of the first ratio arises from the strange quark mass m s (lGeV) = 150 =F 50 MeV. 
This result indicates that SU(3) breaking effects might be significant. 

The B — f V form factors, V = K*, p, <p, associated with the semileptonic decays B — > Vlv and with 
the radiative decays B — > V'y, can be analyzed in a similar way. The semileptonic form factors and 
penguin form factors are defined via the matrix elements, 
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(90) 



respectively, where e denotes the polarization of the vector meson V. I simply quote the results in 
as shown in Fig. El LCSR has been also applied to the B — > p7 weak annihilation |141|. I142j . the 
penguin form factor in the B — > r] transition |143j , and the B — ► fiwy width |141j employing the photon 
distribution amplitude. 



26 




Figure 10: LCSR results for B — ► V transition form factors. 
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Figure 11: Recent determinations of f% and fr> in the quenched approximation. 



4.2 Lattice QCD 

The B meson decay constant and transition form factors, defined as hadronic matrix elements in the 
previous subsection, can be calculated directly on the lattice. For recent reviews on the application 
of lattice QCD to exclusive B meson decays, refer to |144[ I145j . Many results have been obtained by 
different groups using different heavy quark methods, each of which has different systematic errors. For 
example, the UKQCD and APE groups used an 0(a)-improved Sheikholeslami-Wohlert (SW) action, 
a being the lattice spacing, which is defined at the scales of the c quark mass. Outcomes are then 
extrapolated to the h quark mass following the evolution governed by HQET. The Fermilab group 
(FNAL) |146[ I147j identified and correctly renormalized nonrelativistic operators in the SW action, 
such that discretization errors reduce from 0(arriQ) to 0(aAq CD ). JLQCD adopted the action derived 
from non-relativistic QCD (NRQCD). 

Recent determinations of fs and fu in the quenched approximation are summarized in Fig. II II [148 
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Figure 12: Recent lattice results for F + (q 2 ) and F (q 2 ) in the quenched approximation with F + (0) 
F (0). Prediction obtained from LCSR |133j are also shown. 



with the references, from top to bottom, [III], [HO], [T5T] . [TH2] . [THHj . [TBI] . [THH] . [THE] . [ToTj . [T5H], 
and [IBJ] for f B , and [M3,[IHI], CEE], HS2|, CSS], EBU, and [J3E], and \M\ for / D . Values of f B 
derived for a given heavy quark action are consistent. The solid bands, representing the average for a 
particular heavy quark action, are in agreement with each other. The values of fs and fr> in Fig. ^TJ 
are also consistent with the LCSR results in [133 . 

A large source of uncertainty comes from the extrapolation from the scales of the c quark mass to 
those of the b quark mass. There are other subtle issues, such as the scaling violation from discretization, 
and the determination of Jb at rest and at non-zero momentum |160| and from the temporal A Q and 
spatial Ak currents in the matrix element (0|/LJ5(P)) = fsP^- More detailed discussion on the above 
topics can be found |148j . 



Table 2: B and D meson decay constants with Nf = 2. 



Group 


Jb (MeV) /d (MeV) ^ 


Collins99 [KIT] 

MILC'OO dsn] 

MILC'Ol (N f = 2 + 1) PH 
CP-PACS'OO(FNAL) [J32] 
CP-PACS'OO(NR) [El] 
JLQCD PES] 


186(5)(25)( + g u ) ~ 1.26 

191(6)(±?|)(^ 1 ) -1.10 215(5)(^ 3 7 )(io) ^1.08 
1.23(3)(11) 

208(10)(29) 1.11(6) 225(14)(40) 1.03(6) 
204(8) (29) (^ 4 ) 1.10(5) 
190(14)(7) ~ 1.14 



Calculations of decay constants with dynamical quarks have been available, whose results are listed 
in Table [2] [148 . It is observed that Jb is larger in the unquenched theory. The difference between the 
quenched and unquenched predictions depends on which type of the valence chiral extrapolation (linear, 
quadratic or logarithmic |164p is used. Note that fjj may in fact be smaller than the reported value 
due to discretization effects. It is also observed from Table [2] that the dynamical effect for D mesons is 
smaller than for B mesons. 

The B — » 7i transition form factors F + (q 2 ) and Fo(q 2 ) have been calculated in lattice QCD recently 
|165| I166( 11671 116& IT69 . and the results are presented in Fig. H2] [144,. The data show general agreement 
among different groups, within the quoted uncertainties. The agreement is especially good for the form 
factor F + (q 2 ). Note that the quenching effects may be significant, particularly for the form factor F (q 2 ) 
|17()j . The lattice results are available only for large q 2 (smail recoil). Since the bulk of the experimental 
data is located at small q 2 , one needs to extend the calculation to this region at currently accessible 
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iorm factor at maximal recoil, -ti(U) = l 2 {v) = 1 (U) denned m Fq. (|9U|). is given by |144j 

T BA '*(0) = 0.24(5) (3) , 
which is smaller than the LCSR ones |14()[ I171j shown in Fig. 



(91) 



4.3 Hard-scattering Picture 

In the PQCD approach hard dynamics is assumed to dominate in the B meson transition form factors. 
Soft contribution, though indeed playing a role, is less important because of suppression from the 
Sudakov mechanism. Unlike QCD sum rules, soft contribution can not be included into the PQCD 
formalism in a consistent way: if there is no hard gluon exchange to provide a large characteristic scale, 
twist expansion does not hold. Therefore, soft contribution can not be estimated using the same meson 
distribution amplitudes resulting from twist expansion. If it has to be added, it must be introduced 
as an independent input, similar to the treatment in the QCDF approach. The values of these inputs 
usually come from QCD sum-rule or lattice calculations, which, in principle, can contain perturbative 
contributions. Then a double counting of the perturbative contribution, which exists already in the 
one-gluon-exchange diagrams, may not be avoidable. 

It has been explained that the internal h quark involved in the hard kernel becomes on-shell as the 
momentum fraction x of the d quark vanishes in Fig. El 75. . The contributions to the B — ► 7r form 
factor F BlT are then logarithmically divergent at twist 2 and linearly divergent at twist 3 in collinear 
factorization theorem. It has been argued that as the end-point region is important, the corresponding 
large double logarithms a s In 2 x need to be organized into a jet function J(x) as a consequence of 
threshold resummation [73]. This jet function, diminishing as x — > 0, 1, modifies the end-point behavior 
of meson distribution amplitudes effectively. In [§Hj the following approximate form has been proposed 
for convenience, 



J(x) 



2 1+2 T(3/2 + c) 



[x(l -x)} 1 



v^r(i + c ) 

0.3 is determined from the best fit to Eq. 



(EH). 



(92) 

The above expression is 



where the parameter c 
normalized to unity. 

Similarly, the inclusion of fcr also regulates the end-point singularity, and large double logarithms 
a s In 2 fcr are then produced from higher-order corrections. These double logarithms should be organized 
to all orders, leading to fcy resummation IH2j . The resultant Sudakov form factor, constructed in 
Sec. 13.31 [24 , controls the magnitude of k\ to be roughly O(Ams) by suppressing the region with 
fey ~ 0(A 2 ). The coupling constant a s (\/Am B )/7r ~ 0.13 is then small enough to justify the PQCD 
evaluation of heavy-to-light form factors • Note that either threshold or resummation smears the 
end-point singularity. However, to suppress the soft contribution sufficiently, both resummations are 
required, such that the reasonable values of the B — > 7r form factors can be obtained. 

The B — ► 7i form factors F + and Fq are written as, 
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with 



h = 16rm 2 B C F r % J d Xl dx 2 J b 1 db 1 b 2 db 2 ^ + (x 1 ,b 1 )[(j^(x 2 ) - ( pl(x 2 )]E(t^)h(x 1 ,x 2 ,b 1 ,b 2 ) , (94) 
f 2 = lQTxm 2 B CF J dx\dx 2 J bidbib 2 db2^+(xi,bi) 
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4>ts{x 2 ){1 + x 2 rj) + 2rv ( ( x 2 )0t(x 2 ) - x 2 ^(x 2 ) 



+2r 7T <jf 7r E(t ( - 2 )h(x 2 ,x l ,b 2 ,b 1 



(95) 
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E(t) = a s (t)e~ SB ^- s ^ 



(96) 



The explicit expressions of the Sudakov exponents Sb and S n are referred to [21] • The ratio r n is defined 
as rv = mo /tub, where the mass scale mo, related to the chiral symmetry breaking, comes from the 
normalization of the two-parton twist-3 distribution amplitudes (jtff. The contributions from (jjPf are 
of leading power j^H] , and need to be taken into account as mentioned in Sec. 13.21 The first (second) 
terms in Eq. (j9*3j) correspond to Fig. EJa) [Fig. GJb)]. 
The hard function is written as 

h(x 1 ,x 2 ,b 1 ,b 2 ) = J(x 2 )K Q (v^I^2?7m B 6i) 

x [0(&! - b 2 )K (y/x^qniBh) I (y/xtfjm B b 3 ) 

+9(b 2 - bi)K (- v /xirym B 6 2 ) Io (V^"t,b&i)] • (97) 

The jet function J(x) suppresses the end-point behavior of the pion distribution amplitudes, especially 
of the twist-3 ones. The hard scales t are defined as 

= max( x /xirym B , 1/h, l/b 2 ) , 
t (2) = max( v /xIr/m B ,l/6 1 ,l/6 2 ) • (98) 

It is obvious that turning off threshold and Ut resummations with a s fixed, Eqs. (JHH) and (|93|) are 
infrared divergent. 

As stated in Sec. 14.21 lattice calculations become more difficult in the large recoil region of the light 
meson. However, this region is the one where PQCD is applicable |2*4"1 1172j . indicating that the PQCD 
and lattice approaches complement each other. In LCSR 130} 1173] dynamics of the B — > 7r form factors 
have been assumed to be dominated by a scale larger than 0(V Am B ). This is the reason the twist 
expansion into Fock states in powers of l/mris applies to the pion bound state. If this assumption is valid, 
PQCD should be also applicable to the B — > 7r form factors. I emphasize that the "soft" contributions 
have different meanings in LCSR and in PQCD. The soft contribution defined in the former has been 
multiplied by the perturbative Sudakov factor in the latter, such that the soft contribution is large in 
the former, but small in the latter. A good explanation has been provided in |174j . The definitions 
of the "hard" contributions are also different, since the twist expansion has been employed for the B 
meson bound state in PQCD, but not in LCSR. Briefly speaking, the contributions of various orders and 
powers have been organized in different ways in LCSR and in PQCD (also different in LFQCD discussed 
below). Hence, the soft dominance concluded in LCSR does not apply to PQCD |175j . and there is no 
conflict between the basic assumptions in the two approaches. For PQCD to be a self-consistent theory, 
it is only necessary to examine the converegnce of subleading contributions. 

For the B meson wave function, the model j2S] 

$+(x, b) = N B x 2 (l - x) 2 exp 

has been adopted in j^H]- The shape parameter ujb ~ 0.4 GeV has been fixed from the fit to the B —>■ ir 
form factors derived from lattice QCD |176j and from LCSR |130j . The normalization constant Nb is 
related to the decay constant fs = 190 MeV through the relation 

I'dxM*) = £*(*, » = o) = ^. (ioo) 

It is easy to find that Eq. (J9*9*|) has a maximum at x ~ A /m#. The models for the pion distribution 
amplitudes can be found in |106j . 

The relative importance of the twist-2 and twist-3 contributions to F + (q 2 ) has been investigated, 
and the results are listed in Table EH It is found that the latter are comparable to the former, consistent 
with the argument that the twist-3 contributions are not power-suppressed. The approximately equal 
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Figure 13: The B — > 7r form factors / + and /o as functions of q 2 (GeV 2 ). PQCD results for ujb = 0.36, 
0.40, and 0.44 GeV are shown in dots. The solid lines correspond to fits to the lattice QCD results |176j 
with errors. The dashed lines come from LCSR [130 . 



q 2 (GeV 2 ) 


0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 


twist 2 
twist 3 


0.120 0.128 0.138 0.148 0.159 0.172 0.188 0.204 0.223 0.243 0.270 
0.177 0.193 0.210 0.230 0.253 0.279 0.308 0.344 0.385 0.432 0.487 


total 


0.297 0.321 0.348 0.378 0.412 0.451 0.496 0.548 0.608 0.675 0.757 



Table 3: Contributions to f+{q 2 ) from the twist-2 and two-parton twist-3 pion distribution amplitudes. 

weights of the twist-2 and higher- twist contributions to F + have been also observed in LCSR |126j . We 
compare the PQCD results of F + (q 2 ) and F (q 2 ) for q 2 = ~ 10 GeV 2 with those derived from lattice 
QCD |176j and from LCSR (130] in Fig. where lattice results have been extrapolated to the small 
q 2 region. Different extrapolation methods cause uncertainty only of about 5% |177j . The agreement 
at large recoil indicates that ujb ~ 0.4 GeV is a good choice. The fast rise of the PQCD curves at slow 
recoil indicates that perturbative calculation becomes unreliable gradually. The range ojb = 0.40 ±0.04 
GeV, corresponding to F + (0) = 0.30 ± 0.04, has been taken as one of the inputs of the PQCD approach 
to two-body nonleptonic B meson decays. 

The same range of ujb has been adopted in the evaluation of the B — > p transition form factors. The 
results, displayed in Fig.HU are also consistent with those from LCSR |178j at small q 2 . It is found that 
the symmetry relation in Eq. ()109j) below holds very well: A\ is larger than V only by 2% in the large 
recoil region, even after considering the pre-asymptotic forms of the p meson distribution amplitudes 
[106 . Taking the fast recoil limit with rj — > 1 and assuming the asymptotic behavior <p v p = <p a pl the above 
form factors are found to obey the symmetry relations [721 H79j , 

V = At, A 2 = A 1 -2r p A . (101) 

Note that the form factors, treated as nonperturbative objects, are not calculated in [72]. Instead, the 
diagrams in Fig. El under an infrared regularization scheme are regarded as perturbative corrections to 
the relations in Eq. (|1U1J1 . For the application of the PQCD approach to the radiative decay B — ► K*j, 
refer to [THO] . 

4.4 Power Corrections 

SCET provides a systematic framework for discussing power corrections to heavy-to-light transitions 
at large recoil. Here I review the heavy quark expansion of the B meson transition form factors in 
SCET. It will be observed that SCET is a powerful tool of deriving the relations among various form 
factors in the large energy limit. There are three independent form factors associated with decays 
into pseudoscalars, and seven independent form factors with decays into vector mesons. In the former 
case except for the form factors F + introduced in Eq. (|82|) . another one Ft is defined via the matrix 
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Figure 14: The B — > p form factors V, A , A\ and A 2 as functions of q 2 . PQCD results are given in 
dots. The solid lines come from light-cone sum rules. 
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At leading power, the number of independent form factors reduces 



mp) q^ 



(102) 



Assume the energy E 



0(mp) of the final state meson, which is related to the momentum transfer q by E — {rn B + mp — 
q 2 )/(2m B ). The universal functions Ap\E,v,n) and Ay(E,v,n,e), defined via the matrix elements, 
(P(P 2 )|£ T W h v \B{Pi)) and {V(P 2 , e)\^TW h v \B(Pi)), respectively, can be decomposed into their most 
general independent Dirac structures allowed by Lorentz invariance and parity [72*1 1179j . 



Ap\E, v, n) 
A ( °\E,v } n } e) 



2E 4t1UP(E) -2E(v e*) 1 5 ^'(E) . 



(0), 



(103) 



That is, only one universal form factor £p(E) is left for the former, and two form factors, ^(E) and 
ffl(E), corresponding to transversely and longitudinally polarized light vector mesons, respectively, for 
the latter. 

It has been argued [SB.: that the small expansion parameter in SCET should be taken as A ~ \JA/E 
for the heavy-to-light transition form factors. The pion emitted in a heavy-to-light decay at large recoil 
carries momentum scaling like p n ~ (A 2 /E, E, A). This pion is produced in a highly asymmetric state, 
composing of a soft quark with momentum p s ~ (A, A, A) and a collinear jet. This jet must have 
an invariant mass squared p 2 = (p n — p s ) 2 ~ EA = X 2 E 2 . Hence, one has the expansion parameter 
A ~ \J A/ E based on the above kinematical consideration. In this case soft fields carrying momenta of 
order A scale like E(\ 2 , A 2 , A 2 ) and are assigned as being usoft as stated in Sec. 13.21 

One decomposes the matrix element of the heavy-to-light current in full QCD as 



(P(P 2 )\qTQ\B(P 1 )) = tr 
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L 4 2 \ L 4 2 J 

using the projectors (1± /b)/2, /h and /h/A on the large and small components of the quark 
fields, labelled by + and — , respectively. The definition for decays into vector mesons is analogous. The 
functions Ap and Ay with k, I = +, — can again be decomposed as 



A { P(E,v,n) = 2Ef P Kl, (E) 



A^ l) (E,v,n,e) = -IE 7s f?\E) - 2E (v ■ e*) 75 /.[ W) (£) . (105) 



Among the 4 + 8 form factors /p and only 3 + 7 are independent due to the equations of motion 
for light and heavy quarks in QCD and translational invariance, 

q»(P\q lfl b\B) = (m b -m q )(P\qb\B) , 
q»(V\qi^b\B) = -{m h + m q ) (V\q lb b\B) , (106) 

with q = P\ — P2. 

With the SCET expansion of the heavy-to-light currents, it is easy to identity the scaling of the 
above form factors: 

ft + \E)=^°\E) (l + OKA)), 

ft~\E) ~\ft + \E) , ft + \E) ~\ft + \E) , ft~\E) ~ X 2 ft + \E) ■ (107) 

At 0(A), the neglect of \E) leaves 3+6 form factors for pseudoscalar and vector mesons, respec- 
tively. The equations of motion in Eq. (jl06j) give two more constraints at this order. Therefore, one 
has 2+5 independent form factors, implying 1+2 form factor relations. 
At 0(A), the three form factor relations are written as [BJJ] 



2 



Rp = ~ ± V Il = , 9 — J l + 0(a s ,X 2 ) 
b T m B (m B + nip) v 

(i-4V.-*i + 4(i + — H 

L q \ y m B (m B + m v )\ 



ni 



Rw 



m B J \ m B J \ m B J m B \ m B 

r s -fi-4)r, ' 

\ mi ) 



A 



= ^(l + 0(a s ,A 2 )) , (108) 
m B v ' 

for the decays into light pseudoscalar, transversely and longitudinally polarized mesons, respectively. 
As q 2 — > 0, the left-hand sides of the above relations vanish exactly. Other form factor relations, which 
receive O(A) corrections, have been also derived in |59j : 

( i -A) v -( i+ ^) ,ji '=°w- ( i -S) r '- ii =5°w- « io9) 

The above result differs from that in [ST] , where the first equation in (|109|) does not receive an O(A) 
correction. 
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Figure 15: Test of the relations in Eq. (jlOSJ) for B —>■ n and B — > p form factors. The dashed line 
corresponds to the form factors predicted from LCSR 130, 140 . The solid line is the O(A ) result 
including the a s corrections [72j . 



To test the form factor relations in Eq. ()108|) . one substitutes the form factors derived from LCSR 
|1301 1140j into the left-hand sides and compare the outcomes with the right-hand sides. The 0(a s ) 
corrections to the heavy quark limit have been included on the right-hand sides. The results are 
displayed in Fig. ^] [59]. The ratios Rp t T,\\ minus the corresponding values in the symmetry limit are 
shown on the left-hand side. The tensor form factors are evaluated at the scale m&. The grey error 
band reflects the theoretical uncertainty from varying the scale of a s from m^jl to 2m^. The difference 
between the dashed and solid curves is an estimate of 0(X 2 ) or 0(a s X) corrections, which are at most 
3% deviations from zero for q 2 up to 7 GeV 2 . Rp,t,\\ divided by their symmetry limits minus 1 are shown 
on the right-hand side of Fig. It is expected that deviations from are of 0(A 2 ) and/or 0(a s ), since 
the q 2 suppression has been divided out. 

Note that the scaling behavior of the quark fields and of the meson states with A is not sufficient 
to determine the scaling behavior of form factors. Considering soft contribution to the form factor, the 
u quark created in the decay of the b quark carries almost all the energy of the light meson, while the 
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described by the eikonal Lagrangian |179| I181j 

£ cik = q n ^(in-D)q n + 0(1/ E q ), (110) 

where q n (x) = e lEq1l ' x /n /nq(x)/4 are the large components of the light quark spinor field. Even 
for light-cone dominated processes this is an atypical configuration (the preferred one having nearly 
equal momenta of the quark and antiquark). For this reason, although the interaction in Eq. (jllUj) is 
spin-symmetric, the symmetry is not realized in the hadronic spectrum, and there exists no relation 
among the soft contributions to the form factors of pseudoscalar and vector mesons. Furthermore, the 
probability that such an asymmetric parton configuration hadronises into a light meson depends on the 
energy of the meson. Hence, the soft contributions to the form factors are energy- dependent functions, 
whose absolute normalization is not known. 

In this respect SCET applied to heavy-to-light decays at large recoil differs from HQET for B — ► £)W 
form factors. In the case of heavy-heavy form factors, spin symmetry relates pseudoscalar and vector 
mesons, and the Isgur-Wise form factor £(v ■ v') is independent of the heavy quark mass. One obtains 
non-trivial form factor relations beyond the leading order in A, because the heavy quark flavor symmetry 
also relates the initial and final hadronic states. 

Recently, an expansion parameter A ~ A/E for heavy-to-light form factors has been proposed [182J, 
which differs from A ~ yAJ~E discussed above [3T|. l5T?j . For A ~ A/E, the external pion, whose 
momentum scales like p n ~ E(l, A 4 , A 2 ), cannot be built up from the combination of a generic usoft 
momentum p s ~ E(X 2 , A 2 , A 2 ) with a generic collinear momentum p c ~ E(l, A 2 , A). It implies that the 
soft mechanism is strongly suppressed in this picture. For A ~ A/E, the pion momentum scales like 
a collinear momentum. In order to make a light meson out of collinear particles and soft particles, 
one has to require the minus component of the total soft momentum, which would scale like EX, 
to be accidentally small, of order EX 2 . However, this implies a phase-space suppression of 0(A/E) as 
explained in jHE] • It has been expected |182j that under the different choices of the expansion parameters, 
the violations of heavy quark symmetry relations between form factors may start at different power of 
A/E. From the view point of the PQCD approach, the spectator on the pion side is as energetic as 
the collinear particles. That is, all the momenta scale like a collinear momentum, and there is no 
phase-space suppression. 



4.5 Radiative Corrections 

Form factors for heavy-to-light transitions are presumably dominated by nonperturbative QCD dynam- 
ics at small momentum transfer and not computable in perturbation theory. Charles et al. have shown 
that certain symmetries apply to this soft contribution, when the momentum of the final light meson 
is large |179j . These symmetries reduce the number of independent form factors from ten to three as 
shown in Sec. 14.41 The corresponding symmetry relations for the form factors are broken by power cor- 
rections discussed above and by radiative corrections. In this subsection I review the evaluation of the 
symmetry-breaking corrections at first order in the strong coupling constant a s , which are dominated 
by short- distance contributions. The formalism adopted below is referred to as the QCDF approach 

m 

In the absence of a hard spectator interaction shown in Fig. IToT a). the light meson is produced in a 
parton configuration, in which the u quark carries all the momentum of the meson, up to an amount 
of 0(A) in the B meson rest frame. The hard part of the vertex correction in Fig. fTBT b) does not 
respect the symmetry relations, but can be accounted for in perturbation theory by multiplicatively 
renormalizing the current [u n r&„] e ff in the effective theory. A hard interaction with the spectator quark 
shown in Figs. flHT c) and fTET d) allows the meson to be formed in a preferred configuration, in which 
the momentum is distributed nearly equally between the two quarks. 

The soft contribution scales like jTH] 




(111) 
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Figure 16: Different contributions to the B — > P(V) transition, (a) Soft contribution (soft interactions 
with the spectator antiquark q' are not drawn), (b) Hard vertex renormalization. (c,d) Hard spectator 
interaction. 

For the hard contribution, both quarks that form the light meson have momenta of 0(rriB), and the 
gluon in Figs. fToT c) and IToT d) has virtuality of order Am^. The resulting scaling behaviour for the 
pseudoscalar meson form factors is 



K%W « o) ~ < V A ^) — • ( 112 ) 

Therefore, the hard spectator interaction is suppressed by one power of a s relative to the soft contri- 
bution. 

The QCDF formula for a heavy-light form factor at large recoil at leading power in 1/tub is then 
written as 

Fi(q 2 ) = Ci £p(E) + <P B ® Ti <\> P . (113) 

The soft form factor £p(E), defined in Eq. (jl()3|) and represented by Fig. ITflTa). obeys the symmetries 
discussed above. The hard-scattering kernel Ti from Figs. fTHT c) and fTHT d) is convoluted with the light- 
cone distribution amplitudes of the B meson and of the light pseudoscalar meson, for which the endpoint 
divergence together with some finite contribution have been absorbed into the leading soft term. The 
coefficient Cj = 1 + 0(a s ) is the hard vertex renormalization from Fig. 116( b). The correction at 0(a s ) 
from the hard vertex renormalization and from the hard spectator interaction have been obtained in 

To absorb the end-point singularities, the factorization scheme has been defined by imposing the 
condition, 

F+=«p, V=™±™ iu «„ (H4) 

m,B vtiy 

exactly to all orders in perturbation theory, similar to the DIS scheme for inclusive processes. Having 
fixed the factorization scheme, all other form factors can be expressed, for example, as 



2E 

?p 

m B 



l + ^(2-2L) 



+ ^AF , (115) 



M + mp 



1 + ^ L4 + 2L 



^AF T , (116) 

47T 
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L = ^— ln^. (117) 

rriB — 2E tub 



Note L -> 1 for E -»• m B /2 (g 2 -> 0). 

The 0(a s ) corrections to the form factors are given by 



with the quantity, 



AF = AFp , AF T = AFp , (118) 



NcrriB J k + J u 



The theoretical uncertainties in the computation of the hard-scattering correction from the moments of 
the meson distribution amplitudes |183j and from the B meson decay constant are all contained AFp. 

Equation ()113|) has been further elucidated in the framework of SCET [S3]. At leading power in l/m b 
and all orders in a s , a B — ► P transition form factor F can be split into factorizable and nonfactorizable 
components, 

F(E) = f(E) + f F (E), 

f(E) = f -^f [dz[dx[dk + T(z,E,no) 
E z Jo Jo Jo 

x J (z, x, k + , E, [Mo, fi)(p P (x, fi)4> + (k + , fi) , (120) 
f F (E) = C^, /i) . (121) 

Compared to Eq. ()113|) . the hard-scattering kernel has been further factorized into a function T char- 
acterized by the scale m b and a jet function J characterized by the scale /i — y/m^h.. Hence, the hard 
coefficients Ck and T are calculated in an expansion in a s (m&). The jet function J is calculable in terms 
of a matrix element involving a s (y/ Arrib). That is, the contributions characterized by mj and y/m b A 
have been clearly separated. Endpoint singularities arise only in the soft, nonperturbative form factors 
£i(E,fi). The convolution integrals in the factorizable terms are infrared finite in collinear factorization 
theorem. 

I explain the difference between the QCDF formulas based on collinear factorization and those in 
the PQCD approach based on kp factorization. In the former the piece with an end-point singularity 
in collinear factorization theorem has been regularized and absorbed into the soft term / NF . In kp 
factorization theorem the end-point singularity is absent, and both / NF and / F can be formulated into 
the factorization formulas. Since / NF remains in the formulas, the form factor symmetry relations at 
large recoil are still respected in the PQCD approach, which are then modified by the less important term 
/ F . This has been shown explicitly in Eq. (jlOlj) . contrary to the criticism in |72j . It is then realized that 
the definition of soft contributions is in fact ambiguous, depending on the theoretical framework that 
is adopted. In QCDF the soft contribution refers to the one with the end-point singularity in collinear 
factorization theorem (plus an arbitrary infrared- finite piece related to a factorization scheme). In 
PQCD it refers to the one from a large (but arbitrary) coupling constant. Therefore, the hard-scattering 
terms in both approaches [in Eq. (|113|) and in Eq. (|95|) ] also collect different contributions. 

To be more specific, I compare the explicit expression for the form factor F + derived based on 

Eq. 



F + (E) = N J dxdk 
a s (// ) 



2E-m B , 2E 



x k + 



m B m b 
<P^x)<P + (k + ) + C(E, n) C(E, n) , (122) 



with the constant, 



N ° ~ ' (123) 
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Figure 17: Typical 0(a s 
photon vertices. 



contribution to the hard-scattering kernel T/. The crosses represent the 




■ — X s,d 



Figure 18: Typical 0(a s ) contribution to the hard-scattering kernel T t n . 



to Eq. (|93p. The Wilson coefficients satisfy C a = C& = 1 at the tree level. Removing all the Sudakov 
factors and dropping the twist-3 contributions, F + in Eq. (|9*3|) from the PQCD approach reduces to 



F + (E) 



N a s mB / dx\d%2 



V 



X\X 2 X\X 2 



(124) 



Using the variable change rj = 2E /ms — 1, it is easy to identify the three terms in Eq. ()124|) as the 
three terms in Eq. ()122|) in sequence. The third term in Eq. ()124|) with the end-point singularity comes 
from the term 1 in the coefficient 1 + x 2 7] of <j) n in Eq. ()95|) . This piece obeys the large-energy symmetry 
mentioned above. The term x 2 rj, whose x 2 cancels a power of x 2 in the denominator, corresponds to 
the hard-scattering piece in Eq. (|122|) . 

A study of the relative importance of soft and hard dynamics has been done in the framework of 
QCD sum rules |174j . The soft contribution without Sudakov suppression was estimated to be 0.22 
(corresponding to Jb ~ 130 MeV). The soft contribution to fsF Bn obtained in |126j is consistent with 
the above value. It was then shown that the Sudakov effect decreases the soft contribution by a factor 
0.4-0.7, depending on infrared cutoffs for loop corrections to the weak decay vertex. Therefore, the 
soft contribution turns out to be about 0.09-0.15, and smaller than the perturbative contribution about 
0.19. It is then possible that the B —>■ n form factors receive significant perturbative contributions, in 
spite of the large theoretical uncertainty in sum rules, for example, from the variation of the Borel mass 



The QCDF formalism has been applied to the study of the forward-backward asymmetry in the 
rare decay B — > V£ + £~ , where V is a vector meson (184] . Below I discuss the simpler modes B — > V'-f 
The hadronic matrix elements are written as, 



(V-y{e)\Oi\B) = F BV Ti + I d^dvT- 1 (£,v)4> B (£,)(j) V (v 



i , (125) 

where e is the photon polarization vector and the operators come from the effective weak Hamiltonian. 
The soft form factor F BV for the B — >• V transition obeys the symmetry relations in the large energy 
limit. In QCDF the next-to-leading-order hard corrections to the weak decay vertex in Fig. 1171 contribute 
to T/ |187j . These contributions are dominated by scales of 0{m h ) and infrared finite. T/ 7 involves the 
hard scattering of the spectator shown in Fig. IT51 

For the B Vj decay, both the type I and type II contributions can be expressed in terms of the 
matrix element (O7): 

G p 



A(B -> Vi) 



V2 



X^a u 7 + X^a c 7 (V 7 \0 7 \B) 



(126) 



38 




H [GeV] p [GeV] 



Figure 19: Dependence of the branching ratios B(B° — > ^*°7) and B(B~ — > p~7) on the renormal- 
ization scale p, where the leading-order contributions (dotted lines), the next-to- leading-order contri- 
butions including only type-I corrections (dash-dotted lines), and the complete next-to- leading-order 
contributions (solid lines) are explicitly shown. 



where are the products of the CKM matrix elements, the coefficients a? consist of the Wilson 
coefficient C-j and the contributions from the type-I and type-II hard-scattering corrections. It has been 
observed that the leading-order value is enhanced by the T^-type correction. The net enhancement 
of a-j at the next-to- leading order increases the branching ratios as illustrated in Fig. 188 • The 
residual scale dependence for B{B — > K* 0r y) and B(B~ — > p~7) at leading and next-to-leading orders 
is also exhibited. 

The branching ratios of the exclusive radiative decays have been measured to be B(B°—>- K* 0/ ~f) = 
(4.44±0.35) x 1CT 5 and B(B + ^ K*+j) = (3.82±0.47) x 1CT 5 [HH]. For the B -> p 7 decay, only upper 
bound exists. The leading-order results from QCDF have more or less saturated the experimental data, 
and the inclusion of the next-to-leading-order contributions overshoot the data. Note that the transition 
form factor F BV and the distribution amplitudes <$b and 0y are both the independent inputs in QCDF. 
The large predicted branching ratio about B(B° — > K* 0r y) = 7 x 1CT 5 could indicate a double counting 
of hard contributions between the two terms in Eq. ()125|) . Therefore, the basic assumption of QCDF, 
in which the transition form factor is a completely soft object, requires a more careful examination. 

4.6 Light-Front QCD 

In the non-relativistic quark model, wave functions best resemble meson states in the rest frame, or 
where the meson velocities are small. Therefore, the form factors calculated in this model are reliable 
only at small recoil. At large recoil, relativistic effects must be taken into account. A consistent and 
fully relativistic treatment |19Uj of quark spins and the center-of-mass motion can be carried out in 
LFQCD |1911 11921 1193j . This method has several advantages: the light-front (LF) wave function is 
manifestly Lorentz invariant in terms of the momentum fraction variables (in "+" components), which 
is in analogy with parton distributions in the infinite momentum frame. Hadron spin can be correctly 
constructed using the Melosh rotation. The kinematic subgroup of the LF formalism has the maximal 
number of interaction-free generators, including the boost operator which describes the center-of-mass 
motion of the bound state. 

LFQCD has been applied to the heavy-to-heavy and heavy-to-light transition form factors [194, 
195, 11961 1197j . However, the form factors were calculated only for q 2 < 0, whereas physical decays 
occur in the time-like region with < q 2 < (m 8 — rrif) 2 , rriij being the initial and final meson masses. 
Hence, extra assumptions are needed to extrapolate the form factors from the space-like region to the 
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Figure 20: (a) The Feynman triangle diagram, (b) corresponds to the LF nonvalence configuration 
and diagram (c) to the valence one. Filled and empty circles indicate vertex functions and LF wave 
functions, respectively. 



time-like region |198t I199j . Recently, the P — > P transition form factors were calculated in the entire 
range of q 2 [2001120 lj . such that the additional extrapolation is no longer required. This is based on the 
observation |202j that in the frame where the momentum transfer is purely longitudinal, i.e., qx = 0, 
the invariant q 2 = q + q~ covers the entire range of momentum transfer. The price to pay is that, besides 
the conventional valence-quark contribution, one must also consider the nonvalence configuration (or 
the so-called Z graph) in order to maintain covariance. The nonvalence contribution vanishes at q + = 0, 
but is expected to be more important for heavy-to-light transitions near zero recoil |194t 1198] I202[ I203J . 
Prescriptions for treating this nonvalence configuration have been proposed. For example, the authors 
of |201j considered the effective higher Fock state, and calculated the contribution in chiral perturbation 
theory. For a relevant discussion of covariance in the LFQCD framework, refer to [204 . 
Below I mention the other two prescriptions 205} 120b] . Start with the matrix element, 



where q = P\ 



(p'(p 2 )\QYQ\P(Pi)) = F + (q 2 )(Pi + P2T + F-(q 2 ) 9" , 

P2 is the momentum transfer. The form factor F is related to F±(q 2 ) by 

,2 



F (q 2 ) = F + (q 2 ) + 



q 



vrir, — mi 



-F-(q 2 



(127) 



:i28i 



Lp IILpf 

Assume a vertex function Ap |194| 1195] . related to bound state Qq of the meson P. The quark-meson 
diagram in Fig. l2~07 a) gives 



(P'(P2)\Q'Y l Q\P(Pi)) = -j tSi a ^ a ^ 



(2*) 



75 



h + m 3 ) 



m\ + ie 



75 



+ rn 2 ) 



-r 



+ mi) 



p 2 — m\ + ie 



(129) 



wij + ie 

with P2 = pi — q and p 3 = p\ — P\. Consider the poles in denominators and perform the integration 
over the "energy" p{ in Eq. ()129j) . One derives 



(P'(P 2 )\Q'YQ\P(Pi)) 



+ 



[d 3 Pi] 



[d 3 Pi] 




(130) 



s 3 =o 



with the definitions, 



[d 3 pi] = dptd 2 p 1T /(Q47T 3 l[pi 



p = tr[7 5 (^ 3 + m 3 )7 5 (^ 2 + "^h^i + ^1)] 
Si = Pi ~ Pi on , i — 1, 2, 3 , 
Pi (3) = Aon — P3(l)on ' 

PiOTL 



[m 



+ p 2 t)I^p! 



:i3i^ 
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structures Ap are necessary m order to smear the divergences m the integrals. 
Associate 



P1P1P2 



2VP1011 • P2on + mim 2 



R l,2 

iff 3 = V (132) 

2 VP2on ■ P3on ~ m 2 m 3 

with the valence and nonvalence configurations of the Melosh transformation |2()7j Rff\ 2 , which creats 
a state of definite spin (S, S z ) out of LF helicity (Ai, A2) eigenstates |2()1|I2()8[I2()9| . Both of them come 
from the internal structure Ap. One then further makes the substitution, 



A, 



Si + S 3 



Ap/ 

► ^1,3 > Q T Q 

Si=0 ^ 2 + ^ 3 



— > i?» 3 0p, . (133) 

5i=0 



The wave funciton v w , normalized to unity, describes the momentum distribution of the constituents 
in the bound state. 

At last, Eq. (|13U|) becomes 



(P'(P2)\Q'l + Q\P(Pi)) = 2P+H(r) , 

H( r ) = / 0/0 vi \ / 4>p(x,k T )(t)pi(x \kq 



A 2 + k 2 T JA' 2 + k 2 



+ / dx 4> P (x, k T )4> P ,(x', hr) 1 = 1 1 o , (134) 

I A 2 + k 2 T J A' 2 + k 2 J 



with the ratio r = P 2 + /P 1 + and the variables, 

.4 = mix + m 3 (l — x) , A' = m 2 x' + m 3 (l — x') , (135) 

x (x' = x/r) being the momentum fraction carried by the spectator antiquark in the initial (final) 
state in the first term of Eq. ()134j) . However, x 1 > 1 in the second term of Eq. ()134j) indicates that 
the momentum p% of the spectator quark is larger than P' + of the final meson. Therefore, the wave 
function (f> p , plays the role of a fragmentation function in inclusive QCD processes. 
The form factors are then given, in terms of H(r), by 

/ 9n (lTr-)H(r + ) - (lTr+)H(r) , . 

F ± (q 2 ) = ; V +; V - +; V ; , 136 

r + — r_ 



' 2 1 ™2 „2 1 o — rsl „2\ 



r± = - — 7r\m P + m P , — q ±2m P Q(q 
2m 



Op 



where the ratios, 



Q(q 2 ) = ^(m 2 P + m 2 pl -q 2 ) 2 -Am 2 P m 2 P ,/2m P , (137) 

are the solutions of q 2 = (1 — r)(m 2 p — m 2 p ,/r). 

Assume some model wave functions, whose parameters can be fixed from the quark masses, decay 
constants, and other experimental data |2()5j . The numerical results of the B — > 71 form factor F + are 
ploted in Fig. |^ The form factor values are consistent with those obtained in the q + = frame followed 
by an anatylic continuation to the time-like region |211j . It is found that the nonvalence contribution 
to heavy-to-light transitions is negligible in the whole region of q 2 except near zero recoil (q 2 ~ Qmax)- 
In addition, for the same final meson, the nonvalence contributions are smaller when the inital mesons 
are heavier. This conclusion is consistent with those drawn in |194[ 11981 12U2[ !2U3j . 
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Figure 21: The B — > 7r form factor F + compared with the lattice QCD results 
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Figure 22: Non- wave-function vertex(black blob) linked to an ordinary LF wave function(white blob). 



Another treatment of the nonvalence state is to adopt the Schwinger-Dyson equation, which connects 
the embedded state (the black blob in Fig. l2*2*J) to the ordinary light-cone wave function (white blob in 
Fig. |2*2*]) |206j . This connection from one-body to three-body sector can be achieved by introducing an 
operator /C, which in general depends on the involved momenta. It is easy to see that the following 
link between the non-wave-function vertex (black blob) and the ordinary LF wave function (white blob) 
naturally arises, 



(M 2 - M*)V(xi,k Ti ) = J [dyltflTWxiMilV^TjMyjATj) , 



;i38) 



(ml+k 2 ^ 



where M is the mass of outgoing meson and Mq 2 = (mf+kf^) / X\ 
due to the kinematics of the non-wave-function vertex. Note that Eq. 
form as the LF bound-state equation. 

Next step is to remove the four-body energy denomenator using the identity, 



)/(— X2) withxi = 1— £2 > 1 
essentially takes the same 



1 



1 



1 



DaD% D&Dn 



D 9 2 D% 



(139) 



One then obtains the amplitude identical to the nonvalence contribution in terms of ordinary LF wave 
functions of gauge boson and hadron (white blob) as shown in Fig. l2*37 d). Hence, the valence and 
nonvalence contributions can be calculated by menas of the ordinary LF wave functions with the latter 
involving an unknown operator /C. It has been argued that the right-hand side of Eq. ()138|) can be 
approximated as a constant for heavy meson decays in the region with small momentum transfer 206J . 
In contact interaction case, it has been verified that the prescription of a constant operator in Fig. l2*37 d) 
is an exact solution of Fig. l2*3T a). The above formalism has been applied to the D — > K transition form 
factors HI]. 
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5 Two-Body Nonleptonic Decays 

In this section I review progress on the understanding of QCD dynamics in two-body charmless non- 
leptonic B meson decays. Topics related to charmed decays will be arranged in the next section. 
Intuitively, decay products from the heavy b quark move fast without further interaction between them. 
This simple picture, supported by the color-transparency argument [HI], leads to the naive factoriza- 
tion. Although the factorization assumption (FA) |212j gives predictions in relatively good agreement 
with data (apart from the color-suppressed modes), it provides no insight into dynamics. Moreover, 
FA suffers serious theoretical drawbacks as explained below. To improve FA, several frameworks, based 
on different assumptions of the dominant dynamics in exclusive B meson decays, have been developed. 
Amomg these, I will discuss the QCDF, PQCD and LCSR approaches. 



5.1 Factorization Assumption 

To explain the idea of FA, I take the decay B° — > D + tt~ as an example. The relevant effective weak 
Hamiltonian is given by 



V2 



ftcS = —7= V ch V* d 



C^O^ + CMOM , (140) 

with the four-fermion operatos, 

Ox = (db) v _ A (cu) v _ A , 2 = (cb) v _ A (du) v _ A , (141) 

and the definition (qiq2) v±A = 917^(1 ± 75)92- To ensure the renormalization-scale and -scheme inde- 
pendences of physical amplitudes, the matrix elements of four-fermion operators have to be evaluated 
in the same renormalization scheme as that for Wilson coefficients and renormalized at the same scale 
\i. Under FA, the matrix element (0(/x)) is factorized into the product of two matrix elements of single 
currents, governed by decay constants and form factors. The naive factorization was first proved in 
the framework of large energy effective theory |181j . and justified in the large N c limit |213j . For nice 
reviews, refer to |214j . 

In spite of its simplicity, FA encounters three principal difficulties. First, the hadronic matrix element 
under FA is independent of the renormalization scale /i, as the vector or axial- vector current is partially 
conserved. Consequently, the amplitude Ci(/i)(0)f act is not truly physical as the scale dependence 
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contributions have been neglected, lhis may be the reason rA tails to describe the color-suppressed 
decays, such as B — ► J/ipK^*\ Strong phases are essential for predicting CP asymmetries in exclusive 
B meson decays. These phases, arising from the Bander- Silverman- Soni (BSS) mechanism |215j . are 
ambiguous in FA. In this mechanism the c quark loop contributes an imaginary piece proportional to, 

C 2 {t)a s (t) J duu(l - u)6{q 2 u{\ -u)- m\) , (142) 

where q 2 is the invariant mass of the gluon emitted from the penguin. Since q 2 is not clearly defined in 
FA, one can not obtain definite information of the strong phase from Eq. ()142j) . 

The scale problem in FA can be circumvented in two prescriptions. First, one incorporates nonfac- 
torizable effects into the effective coefficients [21^ 12171 I218J: 



af = C^ + dOu) 
Ci(/i) + C 2 (/i) 



,cff 



(143) 



where nonfactorizable terms are characterized by the parameters x%- The /i dependence of Wilson 
coefficients is assumed to be exactly compensated by that of Xi{p) |^14j . However, the renormalized 
four-fermion operator by itself still depends on /x, though the scale dependence of (0(/i)) is lost in FA. 
To next-to-leading order, the Wilson coefficients depend on the choice of the renormalization scheme, 
and it is not clear if Xiif^) can restore the scheme independence of the matrix element. 

In the second prescription, (0(/i)) is related to the tree-level hadronic matrix element via the relation 
(O(fi)) = g(fj,){0)tree- The factor g(/x) is obtained by calculating loop corrections to the weak decay 
vertices. Then schematically one writes 

C(n)(0(fi)) = C(fi)g(fi)(0) tTCC = C cS (0) tiee . (144) 

FA is applied afterwards to the hadronic matrix element of the operator O at the tree level. Since the 
tree-level matrix element (0) tree is renormalization scheme and scale independent, so are the effective 
Wilson coefficients Cf and the effective parameters af expressed by |219l 220J 

af = Cf + Cf(±-+ X i 

af = Cf + Cf^+x*) ■ (145) 

Unfortunately, the extraction of g(fi) from the matrix element is infrared divergent. The divergences 
are usually regularized by considering off-shell momenta p for the external quark lines with p 2 < 0. What 
one has achieved is actually 

C dI = C(ji)g(ji,-p i ,\), (146) 

with p 2 being the infrared cutoff, and A a gauge parameter. Obviously, c eff is subject to the ambiguities of 
the infrared cutoff and the gauge dependence. As stressed in |221j , the gauge and infrared dependences 
always appear as long as the matrix elements are evaluated between quark states. The reason has been 
implicitly pointed out in 222J that "off-shell renormalized vertices of gauge-invariant operators are in 
general gauge dependent". Also, the nonfactorizable contributions are included by introducing more 
free parameters as shown in Eq. (|145j) . These parameters, being process-dependent, then make FA even 
less predictive. The difficulty in predicting strong phases in FA also remains. 

5.2 Q CD Factorization 

An important step towards a rigorous framework for two-body nonleptonic B meson decays in the heavy 
quark limit has been made j2Sl EH I223j . The infrared divergences appearing in the loop corrections 
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Figure 24: Graphic representation of QCD factorization 

to the weak deacy vertices are absorbed into a B meson transition form factor, such that g(fi) can be 
evaluated in terms of on-shell external quarks. In this way, the infrared divergences are regularized 
without breaking the gauge invariance. The nonfactorizable contribution is calculated in the framework 
of collinear factorization theorem, since it is not suffered by the end-point singularity at least at leading 
twist due to the color-transparency argument. The gluon invariant mass q 2 in the BSS mechanism 
can be clearly defined and related to parton momentum fractions in collinear factorization theorem. 
Therefore, the theoretical difficulties in FA are resolved in principle in this QCDF approach. 

The resulting factorization formula for the decay amplitudes incorporates elements both of FA 
sketched above and of the hard-scattering picture. Consider a decay B — > M1M2, where M\ picks up 
the spectator quark. If Mi is either a heavy (D) or a light (it, K) meson, and M 2 a light (n, K) meson, 
QCDF gives the following structure, 

A(B -> M 1 M 2 ) = ["naive factorization"] x [1 + 0(a s ) + 0(A Q c D /m fc )] . (147) 

The 0(a s ) term contains the finite piece of the loop corrections to the decay vertices and the nonfactor- 
izable contributions. The 0(Aqcd/^_b) term collects the power corrections to FA, such as those from 
the annihilation topology. Both terms are supposed to be calculable in a systematic way in QCDF. 
Equation (|147j) . without the chirally enhanced twist-3 contributions, has been proved to all orders in 
a s in the framework of SCET |224| . 

The leading term in l/m& is then written, according to Eq. (j!47|) . as 

(M 1 M 2 \O i \B) = F BMl (0) [ duT^u^M.iu) 

Jo 

+ J d^dudvT n ^,u,v) MOMv)Mu) , (148) 

which is graphically described in Fig. |2U F BMl is a nonperturbative form factor, and 0M; and 0s the 
light-cone distribution amplitudes. The hard kernel T 1 absorbs the finite part of the loop corrections 
to the decay vertices. The hard kernel T n corresponds to the nonfactorizable contributions. It is easy 
to find that Eq. (jl48j) is similar to the leading-power QCDF expression for a B meson transition form 
factor in Eq. (ITT3I) . 

For QCDF, the universal nonperturbative inputs include not only hadron distribution amplitudes, 
but also B meson transition form factors. It has been found that the end-point singularities exist in the 
twist-3 factorization formula for the nonfactorizable amplitudes, and in power-suppressed annihilation 
amplitudes. Because of these end-point singularities, the 0(a s ) and 0(Aq CD /ms) terms in Eq. (|147|) 
turn out to be uncalculable, and their contributions have been parametrized into 

hx^(l + p B e a *) , ln^(l+p^), (149) 

respectively. In fact, the singularities signal the breakdown of factorization. QCDF then contains 
the non-universal and uncontrollable parameters Ph,a and Sh,a- These arbitrary parameters can be 
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Figure 25: 95% (solid), 90% (dashed) and 68% (short-dashed) confidence level contours in the (p, fj) 
plane obtained from a global fit to the CP averaged B — > tttt, Ktx branching ratios. The darker dot shows 
the overall best fit, whereas the lighter dot indicates the best fit for the default hadronic parameter 
set. The light-shaded region indicates the region preferred by the standard global fit jHS], including the 
direct measurement of sin(20i). 



Decay Mode 


Exp. Average 


Default fit 


Fit2 


B° -> tt+tt- 


5.15 ±0.61 


5.12 


5.24 


B ± -> T^TT 


4.88 ± 1.06 


5.00 


4.57 


B° -> 7T 7T 




0.78 


0.94 


B° -» K ± 7T T 


18.56 ± 1.08 


17.99 


18.47 




11.49 ± 1.26 


12.07 


11.83 


B ± 


17.93 ± 1.70 


15.65 


17.88 


B° -> ATV 


8.82 ±2.20 


5.55 


6.87 



Table 4: CP-averaged B — > ixix^Kix branching ratios (in units of 10~ 6 ): data vs. results from the fits. 
The default fit to (p,fj) (returning |Kd>Ac&| = 0.085, 03 = 116° with \ 2 — 4.5) refers to the default 
theory parameter set used in |223j . "Fit2" (returning IK^/V^I = 0.079, 3 = 97°, \ 2 — 1-0) refers to a 
fit without annihilation contributions and chirally enhanced spectator corrections. 



determined, together with the unitarity angles, from the best fit to experimental data |223| P25J. To 
make predictions, QCDF usually presents large theoretical uncertainty due to the arbitrary parameters. 
Another source of theoretical uncertainty comes from the scalar currents, which are proportional to the 
chiral symmetry breaking scale m = m^/(m u + m d ) [23] . 

It is possible to extract the unitarity angle 03 from the data of the CP-averaged B — > ttti and 
B — ► Kit branching ratios. The global fit of the Wolfenstein parameters (p, fj) to the six measured 
B — > mi, Kit world-averaged branching ratios is displayed in Fig. [2U t 227 . The best fits with theory 
parameters in the allowed ranges |223j have x 2 ~ 0.5. The ranges of the strange quark mass and of 
the B meson decay constant are [75, 125] MeV and [170, 230] MeV, respectively. Since a wide range of 
03 is allowed, the result is consistent with the standard fit based on meson mixing and \ V u b\. Figure 1231 
shows a preference for 03 slightly greater than 90°. This result is similar to that from a fit based on 
naive FA |22E| . 

The experimental data and QCDF fits are presented in Tabled |227j . The last two columns come 
from the fitted branching ratios for the default theory parameter set in |223j and the central values of the 
above ranges for m s and and for a second set, where all annihilation effects and chirally enhanced 
spectator interactions are removed. The second set also leads to a good fit without these uncertain 
power-suppressed effects. The normalization of the B — * Ktt modes are sensitive to weak annihilation 
and to the strange quark mass through the scalar penguin amplitude. If 03 is assumed to take values 
around 55° as favored by indirect constraints, the agreement becomes worse for the branching ratios 



46 



-0.75-0.5-0.25 0.25 0.5 0.75 1 

'-'TTTT 



Figure 26: Predicted correlation between the mixing- induced and direct CP asymmetries in the B& 
7t + tt~ decay. See text for explanation of the different curves. 



with significant tree and penguin interference. Note that the inclusion of the annihilation contribution 
weakens the constraint from nonleptonic B meson decays on a global fit of 03 [3Sl [228 . 

The direct and mixing-induced CP asymmetries in the B® — ► 7r + 7r~ decay are also the important 
quantities for extracting the unitarity angles. The time-dependent asymmetry is defined as 

^ T(B°(t)^n+n-)-T(B°(t)^n+n-) 



T(B°(t) -> 7T + 7T-) + T(B°(t) -> 7T + 7T-) ' 

S nn sin(Amt) — C nn cos(Amt) , (150) 



where the direct and mixing- induced asymmetries, 



^irn — 1 i |\ \o. ' 7r7r — 1 i l\ 12 ' V 101 / 



1 - 




2 


1 + 




2 



1 + |A 

respectively, satisfy the relation C% n + S% n < 1. The factor A T7r is given by 

1 + i? c e i<5 e i<t>3 



\ - I \ I p 2j(</. 2 +A0 2 ) _ 2i<p 

' '7T7T ' '7T7T ^ ^ 



1 + # c e i<5 e"^3 



(152) 



with the penguin-over-tree ratio R c — \P/T\ and the strong phase difference between penguin and tree 
amplitudes, 5 = 5p — 5t- We have S nn = sin(20 2 ), if the penguin amplitude is zero. 

The predicted correlation between S nn and C n7T from QCDF is shown in Fig. |2H 227 , where the 
BB mixing phase has been fixed to sin(20i) = 0.78. Each closed curve is generated by specifying the 
theory input and varying 03 from to 360°. The central (dark) curve refers to the calculation of P/T 
with the default theory parameter set, the two neighboring (lighter) curves refer to P/T plus- minus its 
theoretical error without the error from weak annihilation (but including the one from and the 

final (lightest) curves also include the error from weak annihilation. The black part on each curve marks 
the point 3 = 60°; the fat line segment marks the range [40°, 80°] favored by the standard unitarity 
triangle fit with larger 3 to the right of the black part. Note that the Belle data |229j are close to the 
boundary of the physical region. 
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B° 
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6.3; 



8.6; 



7T 



. if* 



0.7 



4.1 
1.5 

+0.6+3.0 
0.5-0.6 



9.2l|? 



< 35 

< 13 
26.5 ±6.1 
16.4 ±3.0 

< 31 

< 3.6 



Table 5: Predictions for the CP-averaged branching ratios (in units of 10~ 6 ), assuming <ft 3 = 70°, 
\Vcb\ = 0.041 and IV^/V^I = 0.09. The first error is due to parameter variations, and the second one 
shows the estimate of the uncertainty due to weak annihilation. The column labeled "default" refers to 
m s = 100 MeV and F 2 = 0. 



The recent BaBar measurement [230J with 90% CL. interval, taking into account the systematic 
errors, are 

= 0.02 ±0.34 ±0.05 , [-0.54, ±0.58], 

= -0.30 ± 0.25 ± 0.04 , [-0.72, ±0.12]. (153) 

It is observed that the QCDF predictions prefer a small C nn , which are close to the upper bound of the 
BaBar data. The predictions for SVtt depend on 02- The determination of the CKM matrix elements 
from the measurement of S nn in the QCDF framework has been performed in |231j . 

The QCDF formalism has been extended to the study of B -> VP modes HSl 1233 • 
The B —* <pK branching ratios were predicted to be around 4 x 10 -6 |232t I233j . which seems to be 
smaller than the experimental data (see Table 7). The reason is that the same set of free parameters in 
Eq. ()149|) has been adopted for the B — > PP and VP decay amplitudes. The annihilation contribution 
is then constrained by the B — > Ki\ branching ratios, and can not help to increase the B — > (f)K 
branching ratios. In a global fit performed in [233] these parameters have been assumed to be different 
for the B — > PP and VP modes. Introducing two independent sets of free parameters, the B — > (f)K 
branching ratios can be fit (due to a larger annihilation contribution) without increasing the B —>■ Kir 
ones. However, the data for the B — > K*tt modes were not included into the global fit. It has been 
noticed |235J that once the B — > K*tt modes are included, the confidence level of the best fit drops to 
below 0.1%. The possible large direct CP asymmetry in the B — ► pit decays |236j also deteriorate the 
fit. 

QCDF has been also applied to flavor-singlet B meson decays, such as B — > K^*'r]^> |237j . It is 
difficult to account for the branching ratios of these modes in FA |238j . The scheme for the rj-rj' mixing, 
with a single mixing angle advocated in |239j . was assumed. The contributions from the b — > ccs 
and b — > sgg transitions through the gluon content of singlet mesons were analyzed carefully. Also, a 
singlet annihilation amplitude, where two gluons radiating from the spectator quark form an 77W meson, 
contributes at leading power. The unknown form factors } (0) were parametrized as 

0)=fl f + J?J ^' (154) 

with P = 7] or 7]' and q = u, d. The decay constants fp and fp are defined through the quark currents. 
In |237j Fx = F^ n (0) and F 2 = or 0.1 were adopted. Combining the above effects, the predictions for 
the CP-averaged B — > K(*'{jf?> , 7T°) branching ratios from QCDF are summarized in Table El 
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Figure 27: Factorization of two-body nonleptonic -B meson decays in the PQCD approach. 
5.3 Perturbative QCD 

As stated before, the PQCD approach to two-body nonleptonic B meson decays is based on kx factor- 
ization. Therefore, the theoretical difficulties in FA can also be resolved but in a way different from 
that of QCDF. The infrared divergences appearing in the loop corrections to the weak deacy vertices 
are absorbed into meson wave functions, such that the infrared divergences are regularized without 
breaking the gauge invariance. The factorizable, nonfactorizable and power-suppressed annihilation 
contributions are calculated in the framework of kx factorization theorem without the end-point sin- 
gularities. The arbitrary cutoffs introduced in QCDF [(""""J 1223j are not necessary, and PQCD involves 
only universal and controllable inputs. The gluon invariant mass q 2 in the BSS mechanism can also be 
clearly defined and related to parton momentum fractions. 

The amplitude for the B — > M1M2 decay is then factorized into the convolution of the six-quark 
hard kernel, the Wilson coefficient, the jet function and the Sudakov factor with the bound-state wave 
functions as shown in Fig. |27| [2U EH| 12321 12331 . 

A = (j) B ® H {6) ® J ® S ® <p Ml ® 0m 2 , (155) 



all of which are well-defined and gauge-invariant. J denotes the jet function from threshold resummation 
discussed in Sec. 12.31 and S denotes the Sudakov factor from fc T resummation discussed in Sec. 13.31 
J (S), organizing the double logarithms in the hard kernel (meson wave functions), is hidden in H 
(the three meson states) in Fig. 123 The partition of nonperturbative and perturbative contributions is 
quite arbitrary. Different partitions correspond to different factorization schemes. However, the decay 
amplitude, as the convolution of the above factors, is independent of factorization schemes as it should 
be. 

The six-quark hard kernel consists of the diagrams with at least one hard gluon |244j . The 
complete set of leading-order diagrams for the B — > Kit decays is displayed in Fig. [2H1 Figures l2*HT a) 
and I28f b). referred to as the factorizable emission, correspond to the leading contribution in QCDF 
[the left-hand diagram in Fig. I24j . Figures |2Hfc) and l2""ff d). referred to as the nonfactorizable emission, 
correspond to the next-to- leading-order contribution in QCDF [the right-hand diagram in Fig. 124*] . 
Figures |2Hfe) and l2""ff f). and Figs. |2""Jg) and I2"""f h) are referred to as the factorizable annihilation and 
the nonfactorizable annihilation, respectively. They are explicitly power-suppressed. However, for the 
physical mass mg ~ 5 GeV, the scalar contribution from the penguin operators, proportional to mo/m^, 
is not really negligible. 

The factorization limit of the PQCD approach at large tub, which is not as obvious as in QCDF, 
has been examined |245j . It was found that the factorizable emission amplitude in Figs. !2HT a) and 
l2"""T b) decreases like m^ 3 ^ 2 as displayed in Fig. [2JJa), if the B meson decay constant fs scales like 
fs oc m B . This power-law behavior is consistent with that obtained in |691 ITS] . Define r as the 
ratio of the magnitude of the nonfactorizable emission amplitude [from Figs. f2HT c) and l2*81f d)] over 
the factorizable one. Figure "29T b). exhibiting r as a function of m_B, indicates that the curve actually 
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Figure 28: Leading-order diagrams for the six-quark hard kernel in the PQCD approach. 
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Figure 29: (a) The factorizable emission amplitude as a function of mg in unit of GeV. (b) The ratio 
r of the nonfactorizable emission amplitude over the factorizable one as a function oi tub- 



descends with m# despite of small oscillation. If parametrizing the ratio as 

iNonfact.l 1 
r = — oc 



Fact. ~ ln a {m B /A) ' ^ 
the best fit to the curve gives the power a ~ 1.0 for A ~ 0.4 GeV |245j . This logarithmic decrease has 
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Figure 30: Plot of vs S n7r for various values of 02 with 0i = 25.5°, 0.18 < R c < 0.30 and 
—41° < 5 < —32° in PQCD. The allowed experimental ranges from BaBar measurment within 90% 
CL. was considered. Dark areas are the regions allowed by PQCD for different 02- 

been confirmed up to = 300 GeV. It implies that the PQCD formalism approaches FA logarithmi- 
cally. 

Surprisingly, the behavior of the ratio r with in PQCD is close to that in QCDF. However, 
the reasonings for achieving the same power counting are quite different. In QCDF the f actor izable 
contribution is assumed to be dominated by soft dynamics, and identified as being of 0(a®). The 
nonfactorizable contribution, being calculable, starts from 0(a s ). Because of the soft cancellation at 
£3 ~ 0(A/m#), the nonfactorizable emission amplitude is dominated by the contribution from the 
region with x 3 ~ 0(1). In this region there is no further power suppression, and one has the ratio, 



In PQCD based on kx factorization theorem |9*5^ I96j. both the factorizable and nonfactorizable con- 
tributions, being calculable, start from 0(a s ). However, the Sudakov factor modifies the factorization 
formulas in the way that a pair of nonfactorizable diagrams exhibits a stronger cancellation as 
increases |246j . It turns out that the ratio r also vanishes logarithmically as shown in Eq. (|156j) . 

I then discuss the applications of PQCD to two-body nonleptonic B meson decays. An alternative 
way to determine 02 is to use the time-dependent CP asymmetry in the B®(t) — > 7r + 7r~ decay, which 
provides two constraints from C nn and S nn for three unknown variables R c , 6 and 02 in Eq. (J150)) . If 
one knows R c and 5, 02 can be extracted from the experimental data of C n7T vs S^. Since PQCD 
gives R c = 0.23tg;c)5 an d —41° < 5 < —32°, the allowed range of 02 at present stage has been fixed to 
be 55° < 02 < 100° as shown in Fig. |3U] |247j . Because the strong phase in PQCD is relatively large 
compared to that in QCDF as explained below, a significant direct CP asymmetry C n7r = — (23±7)% was 
predicted, which could be tested by more precise experimental measurement in the near future PQ [248J . 
The central point of the BaBar data in Eq. ()153|) |230j then corresponds to 02 = 78°. Denote A02 as 
the uncertainty of 02 due to the penguin contribution. For the allowed region of 02 = (55 ~ 100)°, one 
obtains A0 2 = (8 ~ 16)°, implying sizable penguin contributions in the — > 7t + 7t~ decay. The main 
uncertainty comes from the value of \V u b\. 

Here I give a simple explanation for the different phenomenological consequences of the CP asymme- 
tries in two-body nonleptonic B meson decays derived from QCDF and from PQCD. According to the 
QCDF power counting rules [231 EH] based on collinear factorization, the factorizable emission diagram 



^qcdf ~ a s (m B ) oc 



1 



(157) 



ln(m B /A QCD ) ' 
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Figure 31: Plot of R K vs 3 with r K = 0.164, 0.201 and 0.238. 

gives the leading contribution of 0(a®), since the B — > tt form factor is not calculable. Because the 
leading contribution is real, the strong phase arises from the factorizable annihilation diagram, being 
of 0(a s mo/rriB), and from the vertex correction to the leading diagram, being of 0(a s ). For m /mB 
slightly smaller than unity, the vertex correction is the leading source of strong phases. In kp factor- 
ization the power counting rules change 246 . The factorizable emission diagram is calculable and of 
0(a s ). The factorizable annihilation diagram has the same power counting as in QCDF. The vertex 
correction becomes of 0(a 2 s ). Therefore, the annihilation diagram contributes the leading strong phase. 
This is the reason the strong phase derived from PQCD and from QCDF could be opposite in sign, 
and the former has a large magnitude. As a consequence of the different power counting rules, QCDF 
prefers a small and positive CP asymmetry C W7r |227j . while PQCD prefers a large and negative C n7T 
[U 1213 1213 12SH- Significant CP asymmetries are also expected in the B -> Kir j2S], B -> KK 
[252 and B — > pK, ujK |253j decays. The last two modes are especially sensititve to the annihilation 
contributions. 

The CP-averaged B — > Kit branching ratios may lead to nontrivial constaints on the angle 03 0111 ■ 
Introduce the observables, 

Br(B° -> K^ir?) r+ 
K ~ Br(B ± KH±) r ' 

A = A CP (B° - K + n-) R K , (158) 

with the tree-over-penguin ratio tk = \T /P \ for the B — > Kit decays, and the strong phase difference 
between the tree and penguin amplitides, 5 = 5t> — 5 pi. One has [247] 

R K = 1 + r 2 K ± \j{r\ cos 2 3 - A\ cot 2 3 , (159) 

with tk = 0.201 ±0.037 from PQCD |2H], and Aq = — 0.110±0.065 from the recent BaBar measurement 
A CP (B° d -> K+n-) = -10.2±5.0±1.6% |2S0] and the present world-averaged value of R K = 1.10±0.15 

ra. 

With the values 5 pi = 157°, 5t> = 1-4° and the negative cos 5 = —0.91 derived in [247 , one constrains 
the allowed range of 03 within 1 a as displayed in Fig. 

For r K = 0.164: exclude 0° < 3 < 6°. 

For r K = 0.201: exclude 0° < 3 < 6° and 35° < 3 < 51°. 
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18.56 ± 1.08 
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11.49 ± 1.26 


9.1t\i 


B ± - 




17.93 ± 1.70 


17.3 ±2.7 


B°- 
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8.82 ±2.20 


8.6 ±0.3 



Table 6: PQCD predictions for the B — > tttt and Kir branching ratios in unit of 10 6 with 3 = 80°, 
Rb = VP 2 + V 2 = 0-38. 



Decay Channel 


CLEO BELLE BABAR 


PQCD 


(j)K° 


5.5^ ±0.6 11.2^ ±0.14 7.7±iH ±0.8 
< 12.3 8.9±|;| ±1.0 8.1^ ±0.8 


10.21^ 

9 6+ 37 
y.o_ 2 .o 




16ijj;| 26.0 ±9.0 


q i +4.y+U.3 
J - i -3. 9-0.2 




< 16 16.2its 15.5 ± 3.8 


lO.Ot^ ±0.0 






o o+i.y+u.6 






2.8ii;g±0.0 



Table 7: PQCD predictions and experimental data for the B — > <pK^ and B — > X*7T branching ratios 
in unit of 10" 6 for 3 = 80°, # 6 = 0.38. 

For = 0.238: exclude 0° < </> 3 < 6° and 24° < 3 < 62°. 

Taking the central value = 0.201, 3 is allowed in the range of 51° < </> 3 < 129°, because of 
the symmetric property between Rk vs cos0 3 . This range is consistent with that from the model- 
independent CKM-fit in The PQCD predictions for the CP-averaged B — > hit and Kit branching 
ratios are listed in Table 6. 

The leading factorizable contributions involve four-quark hard kernels in QCDF, but six-quark hard 
kernels in PQCD. This distinction also implies different characteristic scales in the two approaches: the 
former is characterized by m^, while the latter is characterized by the virtuality of internal particles 
of order ^ Arris ~ 1-5 GeV |251 I249j . It has been known that to accommodate the B —>■ Kir and 
mi data, penguin contributions must be large enough. In QCDF one relies on the chiral enhancement 
by increasng the chiral symmetry breaking scale to a large value m ~ 3-4 GeV |226j . Because of 
the renormalization-group evolution effect of the Wilson coefficients, the lower hard scale leads to the 
dynamical penguin enhancement in PQCD. The dynamical enhancement of penguin contributions in 
the PQCD approach also appears in the study of B — > VP modes |246| I255| I256| l237j . The predictions 
are listed in Table 7. For a vector meson, the mass mo is replaced by the physical mass my ~ 1 GeV, 
and the chiral enhancement does not exist. Therefore, the ways to account for the B — > VP branching 
ratios in PQCD and in QCDF are different. As stated in the previous subsection, the infrared cutoffs in 
Eq. (|149|) for the B — > VP modes have been assumed to differ from those for the B — >• PP modes [234J. 
A larger annihilation contribution can then help to enhance the B — > cf)K branching ratios without 
increasing the B —>■ Ktt ones. 

At last, the B — > Keta^ decays have been analyzed in the PQCD approach |258t 1259] . However, 
the analysis is not yet complete. 

5.4 Light-Cone QCD Sum Rules 

LCSR has been applied to the B — > tttt decays j^U] and the B — > tttt decays [260J recently. The 
fundamental concept, such as the quark-hadron duality, has been briefly explained in Sec. 14.11 Start 
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Figure 32: Sum rule calculation of the B — > txtx decay. The shaded oval region denotes nonperturbative 
input, the pion distribution amplitude. The other pion and the B meson are represented by the currents 
JW(p — k) and j( B \p — q), respectively. The square represents the four-quark operators Oj. 

with the correlation function, 

F u (p,q,k) = J d'xe-^-^J d%e^- k >y(0\T[4l\y)O(0)4 B \x)Mq^ , (160) 

where the external pion state has been specified, the two interpolating currents, 

jj$ = u^ulbd , J ( 5 B) = m b bi~/ 5 d , (161) 
are for the pion and for the B meson, respectively, and the relevant operators are 

Ox = (di»(?zm) , dx = (dr^ u ) ( ErM y 6 ) ■ (162) 

The configuration is illustrated in Fig. with an unphysical momentum k coming out of the weak 
vertex. This momentum was introduced to prevent the B meson four-momenta from being the same 
before (ps = p — q) and after (ps = P) the decay. Then the continuum of light states will not enter 
the dispersion relation of the B meson channel. 

Take p 2 = k 2 = q 2 = 0, and consider the region of large spacelike momenta, 

\{p-kf\ ~ \{p-qf\ ~ \P 2 \ » A^ CD , (163) 

in which the correlation function is explicitly calculable by means of OPE. The decomposition of the 
correlation function in Eq. ()16U|) in the independent momentum structures contains four invariant 
amplitudes, 

FP = (p - k) v FM + q u F[° ] + KPP + e^ Xp qYk p FP , (164) 

for the operators O = Oi,0\, where only the amplitude F^ is relevant. The procedure to derive a 
double dispersion relation is as follows |261j . One first makes a dispersion relation in a pion channel of 
momentum (p — k) 2 and applies the quark-hadron duality for this channel. Thereafter, to extract the 
physical B meson state, one performs an analytical continuation of the invariant mass P 2 to its positive 
value, P 2 = m 2 B . This procedure is analogous to the one in the transition from spacelike to timelike 
form factors. Finally, a dispersion relation in the B meson channel of momentum (p — q) 2 is derived, 
together with the application of the quark-hadron duality [20]. 

At the diagrammatical level, there are four topologically different contributions to the correlation 
function in Eq. (|16U|) . corresponding to four possible combinations of u and d fields in the pion distribu- 
tion amplitude (0 | u a {z\)dp(z2) \ 7r~), with z\ = or y, and z 2 = x or y, a, (3 being the spinor indices. 
Drawing the quark diagrams, one finds that each contribution yields a B — > tttt matrix element with 
a certain quark topology: emission (zi = 0, z 2 = x), annihilation (zi = 0, z 2 = y), penguin {z,\ = y, 
z 2 = x) and penguin annihilation (z\ = z 2 = y). So far, only the emission topology F^ has been 
calculated [20] . 

For the matrix elements of Oi, the factorizable diagrams are those, in which the quarks of the 
heavy-light currents do not interact with the quarks of the light-quark currents. A typical example is 
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p-q k 

(a) (b) 

Figure 33: Diagrams corresponding (a) to the leading-order of the correlation function in Eq. (|16U|) 
for O = 0\\ (b) to the higher-twist soft-gluon nonfactorizable contribution for O = 0\. Solid, dashed 
and wavy lines represent quarks, gluons, and external momenta, respectively. Thick points denote the 
weak interaction vertices, and ovals the pion distribution amplitudes. The cross represents another 
attachment of the gluon. 

shown in Fig. I3*37 a). To calculate the factorizable contribution, one inserts an intermediate vacuum 
state between the weak currents of the operator 0\. Equation (jl6(J|) is then converted into a product 
of two disconnected two-point correlation functions, 



F%\p,q,k) = \ij Ae*" fc) ' J (0 | T[j%(y)d(0) W u(0)} | 0)J 

x(i[ d A x e iip - qyx (0 | r[«(0)7"6(0)ji B) (ar)] | ir~(q)) J . (165) 



The analysis then reduces to that of the B — ► 7r transition form factors. 

For the operator 0\, nonfactorizable corrections to Eq. ()16())1 . appearing at a two-gluon level, are 
negligible. In the case of Ox, nonfactorizable effects start at the one-gluon level. The relevant correlation 
function Fj)°^ receives contributions of hard gluon exchanges, whose 0(a s ) examples are shown in 
Fig. EU These two- loop diagrams, not yet calculated because of their complexity, are very important: 
they give the scale dependence of the matrix element, which partially compensates the scale dependence 
of the Wilson coefficients in the effective weak Hamiltonian. Moreover, the analytic continuation of 
these two-loop contributions in P 2 generates an imaginary part, which is essential for predicting the CP 
asymmetries in the B — > tiit decays. The diagrams in Figs. l3~3T a) and l3"3T b) correspond to the corrections 
to the weak decay vertex in the literature. Figure l3~3T c) corresponds to the hard spectator contribution. 
However, as explained before, the soft and perturbative contributions in the QCDF, PQCD and LCSR 
approaches all have different definitions. 

There is another type of nonfactorizable effects, which comes from the diagram with on-shell gluons 
being emitted from the quarks of the pion current and absorbed into the pion distribution amplitude, 
as shown in Fig. l33T b). In terms of the light-cone expansion these contributions are of higher twist, 
starting from twist 3. It has been argued that the higher-twist nonfactorizable effects are suppressed 
by a power of 1/tub compared to the twist- 2 factorizable amplitude [20 . To quantify the magnitude of 
the nonfactorizable effect from the twist-3 pion distribution amplitude, the ratio has been introduced, 

\ E (B° -> 7T + 7T-) _ 4 5l) (5° - 7T+7T-) 



Xe was estimated to be 



4 0l) (SS - 7T+7T-) 



\e(B® — > 7r + 7r~) = 0.05 -J- 0.15 GeV , (167) 



indicating that nonfactorizable soft corrections from the twsit-3 pion distribution amplitude appeared 
to be numerically small (~ 1%). 

Recently, a piece of higher-order contribution to the B — > tttt decays from the chromomagnetic 
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(a) 



(b) 




Figure 34: Diagrams corresponding to the 0(a s ) nonfactorizable contributions to the correlation func- 
tion in Eq. ffToTljl for O = 6 X . 



dipole operator 0$ g (gluonic penguin) has been evaluated in LCSR |175j . Similarly, consider the ratio 



o 



7T 1 7T 



r (Os g )(B» _> vr+vr-) = '- , (168) 

4 0l) (^ -^-) 

which determines (up to the known Wilson coefficient Cs g ) the gluonic-penguin correction to the fac- 
torizable B — > tttt decay amplitude. The result 

r {Oss) (B d -> tt+tt") = 0.035 ± 0.015 , (169) 

is of the same order of magnitude as in Eq. ()167|) . Though the impact of gluonic penguins on the 
B® — > 7T + 7r _ mode is very small, it might be noticeable for the B® — > 7r 7r mode |175j . 



6 Charmed Decays 

6.1 Transition Form Factors 

B and D meson hadronic matrix elements play a crucial role in the determination of the CKM matrix 
elements and in overconstraining the unitarity triangle of the Standard Model. The B — > D^*' transitions 
are defined by the matrix elements, 



(L>(P 2 )|&(0) 7m c(0)|S(P 1 )} 
(D*(P 2 ,e)\b(0)^c(0)\B(Pi)) 

(D*(P 2 ,e)\b(0)lAV)\B(Pi)) 

with the velocity transfer 77 = V\ • v 2 , v± 
£a 2 > £a 3) an d £y satisfy the relations in the heavy quark limit, 



y/m B m D [£+(rj) {v x + v 2 ) ll + £-(v)( v i ~ ^2)^ , 

-Usiv)^* ■ Viv 2/J ] , 
i y /m B m D ^ v (r])e flua ^ 'elv 2a v 1/3 , 



(170) 



P\/ m B an d v 2 = P 2 /m D (*). The form factors £+, 

0, (171) 



56 



. 8 
. 6 
. 4 
0.2 



1.25 1.5 1.75 2 2.25 2.5 2.75 3 

Figure 35: Isgur-Wise form factor £(77) for different Borel parameters. 
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Figure 36: 1/m correction form factors pi(r?) (solid curves), P2(v) (dashed curves), and piirj) (dashed- 
dotted curves) in the unit of 100 MeV. 

where £ is the Isgur-Wise (IW) function |262j . 

The 0(l/m c ) corrections introduce four new functions 262], and the 0(l/mb) corrections do not 
[214 . Taking the matrix elements of the vector current as an example, one has 




(172) 



with the mass difference e = rriB — rrib = mp — m c = mo* — fn c , all of which are equal up to 1/m 
corrections. The Luke's theorem |214l I263| I264j leads to the values of the subleading form factors at 
zero recoil, pi(l) = 0, P2(l) = 0. It is essential to examine these 1/m corrections, especially those from 
l/m c . If they are modest, HQET will be self-consistent and useful. 

The form factors in Eq. ()170|) have been calculated using QCD sum rules at finite m^, m c in |265| l266j. 
HQET sum rules for the IW function were derived in |267t 1268] . which coincide with the limit m^^ — > oo 
of the QCD sum-rule calculation. The above results for finite and infinite masses were compared in 
[269 . There are two alternative ways to obtain sum rules for the subleading form factors pi. One can 
either expand the known finite- mass QCD results to the first order in 1/m, or start from the HQET 
Lagrangian and currents in the first power in 1/m. It has been noticed that the slope of £(r?) near rj = 1 
depends on how to model the contribution from the continuum state on the hadron side |267[ 1268] f270j. 

The sum rule for £(77) is well known |2671 1268] . whose results for several Borel parameters are shown 
in Fig. 1351 The predictions for the 1/m-correction form factors pi t 2,4,{v) are presented in Fig. ESI t 270j . 
The curves of pi^iv) are pinned at the origin by the Luke's theorem, and increase for 77 > 1. p^rj) 
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estimate their accuracy. Hence, it is expected that the U(l/m) corrections to the B — > D y > transition 
form factors are at most 10%. The above analysis provides the behavior of the B — ► D transition form 
factors with the velocity transfer, but the values at zero recoil, which are specially important for the 
extraction of the CKM matrix element \V c b\, remain intact. 
The IW function near the zero recoil is parametrized as 

£{V) = F B ^ DW (l)[l - p 2 DM (v - 1) + c DW (r) - l) 2 + 0(( V - l) 3 )] . (173) 

Including the 1/m and radiative corrections, the normalizations Fb->d(1) = 0.98±0.07 and = 
0.91 ± 0.03 have been derived in |271j . These normalizations can also be studied in lattice QCD, and 
the results are 

F B ^ D (1) = 1.058 ± 0.016 ± 0.003+g-^f fI7I\ . 

iW(l) = 0.9130lg:KSS m • (174) 

The errors come from fitting, matching lattice gauge theory and HQET to QCD, lattice spacing depen- 
dence, light quark mass dependence and the quenched approximation. The above values agree with those 
from other methods, such as non-relativistic quark models |274j and a zero-recoil sum rule |275M276] . For 
a QCD sum-rule analysis of the subleading form factors involved in the semileptonic decays B — > D\lv 
and B -> D\lv, refer to |277| . 

The LFQCD formalism has been applied to the B — > D* form factors, whose results can be found 
in j278] . 

The PQCD formalism for B — > transitions has been developed recently [279J, which applies 
under the hierachy, 

m B > m DW > A , (175) 

with m D (,) being the meson mass. The relation m# ^> m D (») justifies perturbative evaluation of 
the B — ► form factors at large recoil and the definition of light-cone meson wave functions. 
The relation m D o) ^> A justifies the power expansion in the parameter A/m D (*). Equation ()175|l . 
corresponding to the heavy quark and large recoil limits, may not be realistic. Nevertheless, an attempt 
to construct a self-consistent theory under this hierachy is worthwhile. 

It has been argued that the wave function for an energetic D^*> meson absorbs collinear dynamics, 
but with the c quark line being eikonalized. That is, its definition is a mixture of those for a B meson 
dominated by soft dynamics and for a pion dominated by collinear dynamics. The behavior of the 
heavy meson wave functions under Eq. ()175j) has been examined. For A/rns, A/m^.) 1, only a 
single B meson wave function (f>+(x) and a single meson wave function (p D (*)(x) are involved in 
the B —> D(*} form factors, x being the momentum fraction associated with the light spectator quark. 
Equations of motion for the relevant nonlocal matrix elements imply that <f>+(x) and <fi D (*)(x) exhibit 
maxima at x ~ A /tub and at x ~ A/m D M, respectively. To proceed a numerical analysis, the simple 
model [2791 . 

0dm(z) = -j^=f D (*)x(l-x)[l + C DM (l - 2x)\ , (176) 

has been adopted. The free shape parameter C D (» is expected to take a value, such that d m has a 
maximum at x ~ A/m D (*). The intrinsic b dependence of the D meson wave function was not included. 

The free parameters C D (*) can be fixed by fitting the leading-power PQCD predictions to the mea- 
sured decay spectra |280t I28l| l2~82 . With the normalization given in Eq. (|174j) . the linear and quadratic 
fits to the data give j2HDl I2HH 

p 2 D = 0.69 ±0.14 , c D = 0, p 2 D * = 0.81 ±0.12 , c D * = , 

& = 0.69j£g, c D = 0.00^, (177) 
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Figure 37: (a) [(b)] £ as a function of the velocity transfer from the B — ► D^*Hu decay. The solid lines 
represent the central values, the dashed (dot-dashed) lines give the bounds from the linear (quadratic) 
fits. The circles correspond to C D (,) = 0.5, 0.7, and 0.9 from bottom to top. 

respectively. Choosing the decay constants fs = 190 MeV and fo = /d* = 240 MeV, it has been 
found that Co ~ Co* = 0.7 ± 0.2 leads to an excellent agreement with the data as exhibited in 
Fig. |37| For these values, the corresponding meson distribution amplitude shows a maximum at 
x ~ 0.36, consistent with the expectation. The rough equality of Co and Co* hints that the heavy 
quark symmetry holds well. It has been shown that the leading PQCD factorization formulas B — > 
transitions indeed respect the heavy quark symmetry. 

6.2 B -> Dtt 

The recent measurement of the B° — > D°tt° branching ratio reveals interesting QCD dynamics. In naive 
(or generalized) FA, nonfactorizable effects are parameterized through the phenomenological coefficients 
a| ff in Eq. (j!45|) . which depend on the color and Dirac structure of the operators, but otherwise are 
postulated to be universal |212|I214| 1283] . Class-1 and class-2 decay topologies refer to the cases, where 
a charged and a neutral final-state meson are produced from the four-quark operators, respectively. For 
instance, the decay B° — > D + n~ is a class-1 process, in which the charged pion is generated at the 
weak vertex, whereas B° — > D°7r° is a class-2 process, in which the D° meson is directly produced. The 
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5.1 ± U.4 ± U.b 


a.Y4Ao:32 


± U.bb 


B° - 


-> D*°7T 


o 7+O.8+O.5 
z - ' -0.7-0.6 


9 9O+0.59 
z - zu -0.52 


±0.79 


B° - 


-> D°r/ 


1.4lg;t ± 0.3 






B°- 


-> D*°r/ 


2.0 +0 ,| ± 0.4 






B°- 


-> D ^ 


1.8±0.5lgj 






B°- 


-> L>* cj 


3.1±i;f d= 0.8 







Table 8: Data (in units of 10" 4 ) of the B° -> £)(*)°M° (M = vr, 77, w) branching ratios. 



corresponding amplitudes are then expressed as 

^#-^°) = l ^V cb V: d (m 2 B -ml)f D Ft(m 2 D )a 2 (Dn), (178) 

where the coefficients have the orders of magnitude a\(Dir) ~ 0(1) and a 2 (D-K) ~ 0(1/N C ). The isospin 
symmetry then implies 

-> DV") = A(5° -> D+7T-) + ^2 -> £>°7r°) . (179) 

Within errors, the class- 1 decays B° — ► D^ + M~ with M = tt, p,ai, D s , D* can be described using 
a universal value |ai| ~ 1.1 ± 0.1, whereas the class-2 decays B — > K^*'M with M = J/ip,ilj(2S) 
suggest a nearly universal value |a 2 | ~ 0.2-0.3 |284j . The wide range of | a-2 1 is due to the uncertainty 
in the B — > form factors. Moreover, the class-3 decays B~ — > D^*'°M~ with M = ir,p, which are 
sensitive to the interference of the two decay topologies, could be explained by a real, positive ratio 
02/01 ~ 0.2-0.3, which seemed to agree with the determinations of |ai| and | a-2 1 from other modes. The 
observed branching ratios of the color-suppressed modes are listed in Table |H] [285, 286J. The parameter 
o 2 extracted from TableHfalls into the range of |o 2 (£>7r)| ~ 0.35-0.60 and \a 2 (D*7r)\ ~ 0.25-0.50 (2BZj - 
The phases of 02/01 are 59° for the Dn system and 63° for D*n [287 , implying sizeable relative strong- 
interaction phases between class- 1 and class-2 B — ► D^'ir decay amplitudes [288, 289J. These results can 
be regarded as a failure of naive FA: the parameters 02 in different types of decays such as B — > n 
and B — ► K^'J/tp differ by almost a factor 2 in magnitude, indicating a strong nonuniversality of 
nonfactorizable effects. 

If the c quark is treated as a massive quark, QCDF does not apply to the class-2 decays B° — > 
_D(*)°M°, because of the uncancelled end-point singularities. Therefore, the magnitude and phase of the 
a 2 (D^M) parameters are not calculable in QCDF. However, these decays are calculable in the PQCD 
approach based on kx factorization theorem, in which the end-point singularity does not exist. From 
the power counting rules proposed in Eq. (|175|) [279 , it has been shown that the relative importance 
of the different topologies of diagrams for the B — > Dtt decays is roughly 

emission : nonfactorizable ~ 1 : — — , (180) 

m B 

which approaches 1 : A /tub as the D meson mass mn reduces to the pion mass of O(A). Since the 
factorizable and nonfactorizable diagrams contribute to the parameters ai and 02 in PQCD, respectively, 
the ratio |a2|/oi ~ 0.5 is obtained. Moreover, the imaginary nonfactorizable amplitudes determine the 
relative phase of the factorizable and nonfactorizable contributions, which is about —57°. 

To obtain the above results, the meson wave function determined from the semileptonic B — > 
D^*nv decay in the previous subsection has been adopted. Therefore, there is no free parameter in the 
above calculation. The PQCD predictions for the B — ► Dn branching ratios [290J, 

B(B' -> D Q n-) ~ 5.5 x 10~ 3 , 
B(B° -> D+tT) ~ 2.8 x 10~ 3 , 

B(B° -> D\°) ~ 2.6 x 10~ 4 , (181) 
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hint large final-state interaction was drawn from the analyses based on rA [287, 288, 292, 293J. Hence, 
it is more reasonable to conclude that experimental data just hint large strong phases, but do not tell 
which mechanism generates these phases |294j . From the viewpoint of PQCD, these strong phases are 
of short distance, and produced from hard kernels. For the application of the PQCD approach to other 
charmed nonleptonic decays, refer to [295J. 



6.3 Decays into Charmonia 

The B meson decays into charmonia, such as B —>■ J/ipK, are theoretically notorious. As defined above, 
this mode belongs to the color-suppressed class-2 decays, in which one expects large nonfactorizable 
contributions. Experimentally, the measurement of the branching ratios 



B(B + -> J/i)K + ) = (10.1 ±0.3 ±0.5) 
B(B° -> J/i/jK°) = (8.3 ± 0.4 ± 0.5) ) 



x 10~ 4 , 

xlO" 4 , (182) 



provides the information of the coefficient a%(J /ipK), being of order 0.20 — 0.30 [284J. However, it is very 
difficult to understand such a large a 2 {J /ipK), and also those for other similar decays into charmonia: 
the observed branching ratio differs from the naive FA prediction by at least a factor of 3. In this 
subsection I discuss the attempts made in the LCSR and QCDF approaches, and explain why they fail. 

The QCDF method is usually not applicable, if the emitted meson is heavy. Take the B° — > D°7r° 
decay as an example. Since the D° meson is not a compact object with small transverse extension, it 
will strongly interact with the (Bit) system, such that the factorization breaks down. The parameter 
a 2 {'KD) has been roughly estimated in [UU] by treating the charmed meson as a light meson, which is 
certainly a dubious approximation. Fortunately, QCDF is applicable to the B — > J/ipK decay, because 
the transverse size of J/ip becomes small in the heavy quark limit. However, a recent study [297 
indicates that the leading-twist (twist-2) contributions from QCDF are too small to explain the data. 

The authors of 298j then calculated twist-3 contributions from the diagrams in Fig. EH The result 
for Figs. [3*8T a)- f3*8T d) is free of the end-point singularity. The contribution from Fig. EEt e ) and l3*8T f) is 
written, up to twist 3, as 

fn = fi? + fi? + ---, (183) 

(2) 

where the superscript denotes the twist of the kaon distribution amplitude. The expressions for f n 
and ffi are 

/}? = A 4 r,J K {\ 2 — T d Xl ^ t dx 2 ^ (X2) i 1 dx 3 , (184) 



N c F+ K (m 2 ^)m 2 B 1 — r Jo ± x\ Jo " x 2 Jo J x 3 
{3) /2m \ 4vr 2 f K f B fi dx x [ l dx 2 , , f 1 dx 3 (j) a K (x 3 ) 

f » = N: F^(m^)m% Jo ~ ^ L J ^6(1^' (185) 

with the mass ratio r = rn 2 ji^/m B and the B — > K transition form factor F+ K (m 2 ^) being evaluated 
at the J /if) meson mass. The twist-2 and one of the two-parton twist-3 distribution amplitudes are 
given, in their asymptotic form, by 

4>k( x ) — 6x(l — x) , 4>k( x ) = — x) , (186) 

(2) 

respectively. It is observed that /)/ is finite, because the potential logarithmic divergence cancels 
between Fig. l3*8T e) and lHBT f). ffj is singular, since only the potential linear divergence cancels, leaving 
the logarithmic one. 

To estimate the twist-3 effect, the divergent integral has been parametrized as 



X= I 1 = In f^±|+p, (187) 



ox 3 V A 



QCD , 
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Figure 38: Feynman diagrams for nonfactorizable corrections to B — > J/ipK. 
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Figure 39: The coefficient \a,2(J/if>K)\ vs. the phase of the parameter r. Solid and dashed curves are 
for \p\ = 6 and 3, respectively. The upper and lower solid curves are for r = rn 2 j^/ml and m^^/m 2 B , 
respectively, and likewise for the dashed curves. 



as in Eq. ([149)1 . with p = \p\exp(i5) being a complex random number. The variation of \a2(J/ipK)\ 
with the arbitrary parameter p is exhibited in Fig. EH1 L 298j, which implies that \ci2(J/ipK)\ can fit the 
data, only when \p\ is almost as large as 6. For such a huge subleading contribution, the self-consistency 
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A worse situation has been observed: QCDr breaks down tor the B — > x&K deacy even at twist 2 
Similarly, the function fjj is obtained by computing the two spectator corrections in Figs. EHe) 
and I38f f). whose contribution is given by 



fn = A l ^rk, K { B \ i ■> I dxj^-^- f dx 3 ^^- [ dx 2 (p Xcl (x 2 ) 



Xcl' 



N c F£ K (ml Jm| 1 -r Jo Xl Jo x 3 Jo Xc 



1 2r 



x 2 x 3 (l 



with the mass ratio r = m Xcl /mB- Obviously, the integral over x% in Eq. (j!88|) gives logarithmic 
divergence. Therefore, QCDF breaks down even at leading twist. This is different from the B — > J/ipK 
decay, which does not have logarithmic divergence at leading twist |297( I298j . The reason is that the 
logarithmic divergences arising from the contribution of the vector and tensor currents are cancelled 
out in the B — > J/ipK decay, while there is no such cancellation for the B — > XciK decay. 

The logarithmically divergent integral has been parametrized in the same way as in Eq. (|187|) . In 
[299 , the parameter X is chosen as X ~ 2.4 to make a rough estimate. For <f) Xcl (x) = 5(x — 1/2), 
B{B — > XciK) = 0.16 x 10~ 4 was obtained |299j . The measured branching ratio |301| 

B(B° -> XciK ) = (5.4 ± 1.4) x 10" 4 , (189) 

is about thirty times larger than the theoretical prediction. Choosing X w 2.4, instead of around 6 
298 , the QCDF prediction for the B — ► J /ipK branching ratio is also too small. 

The end-point singularity becomes more serious in the B + — > XcoK + mode. In the previous cal- 
culations, the contribution of the four vertex diagrams in Fig. EH1 is infrared safe. However, for the 
B — »• XcoK decay [299,, these four diagrams produce infrared divergences. The QCDF predictions for 
the B — > tjPK branching ratios, thought infrared safe, are also too small [300 . The above analyses 
indicate that it is difficult to apply QCDF to B meson decays into charmonia. 

I now turn to the LCSR approach. For the B — > J/ipK decay, the relevant operators are 



O x = (H»(sF*&) , Oi = ( cvj±c J [sV^b ) . (190) 



In the factorization limit, the matrix element of 0\ vanishes, and the matrix element of the operator 
0\ can be factorized into 



(JMp,e)^(g)|C> 1 |5(p + <?)> = (J/^p,e)\cT^c\0)(K(q)\sm\B(p + q)) 

= 2e- ? m w / J/ /f(m^), (191) 



with Jjjfy being the J/ip meson decay constant. The B — > K transition form factor F^ K (m^) = 
0.55 ±0.05 [26 lj was obtained using LCSR in a way the same as for the B —>■ n form factor. Evaluating 
the B — > J/tpK branching ratio with the next-to-leading-order Wilson coefficients, one arrives at [261 

B(B -»■ J/i)K) fact = 3.3 x 10~ 4 . (192) 

This value, representing a prediction from FA, is too small compared to the experimental data in 

Eq. (HH2D- 

Similarly, the nonfactorizable contribution associated with the operator 0\ starts at the one-gluon 
level, and that with 0\ starts at the two-gluon level. One then considers only the matrix element of 
Ox, which is expressed as 

(Jlil>K\d x (ji)\B) = 2e • ?m w / w F + M (/i 2 ) . (193) 
To calculate the factor F^ K , one studies the correlation function, 

F„(p,q,k) =i 2 I d A xe- l ^ x [ d%e^ p - k > y (K(q)\T[4 J ^\y)O(0)4 B \x)]\0) , (194) 
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The results are F+ K ^ 3 \iJ< b ) = 0.003 — 0.0055 from the twist-3 kaon distribution amplitudes and 
F+ K ^\fAb) = 0.006 — 0.0012 from the twist-4 kaon distribution amplitudes with the scale fib — mb/2. 
The wide range for F+ K is attributed to the variation of sum-rule parameters. Combined with the 
factorizable contribution, one derives the effective parameter 02, 

a 2 ~ 0.14-0.17|^, (195) 

which is still too small. 

Since the PQCD approach based on factorization is free from end-point singularities, it is ap- 
plicable to B meson decays into charmonia. The formalism for the color-suppressed nonfactorizable 
amplitude is the same as that for the B — > D°7r° decay in the previous subsection. The B — > J/t^K^ 
decays have been analyzed recently [302, 303], and the predicted branching ratios, together with the 
experimental data, are listed in Table |H1 Briefly speaking, the measured B° — > J/^>K° branching ratio 
|301j is employed to determine the unknown J / ip meson wave function. This wave function is then used 
to predict the B° — > J/^/K*° branching ratio, whose consistency with the data is obvious. 



Mode 


Belle [304J 


BaBar [30 lj 


PQCD 




7.9 ±0.4 ±0.9 


8.3 ±0.4 ±0.5 


8.3 




12.9 ±0.5 ±1.3 


12.4 ±0.5 ±0.9 


13.37 



Table 9: PQCD results of the B — ► J/ipK^ branching ratios in unit of 10 4 with the J/ip meson decay 
constant fjm = 0.405 GeV. 



7 Other Topics 

7.1 Final- State Interaction 

A strong phase can come from short- distance dynamics and from long-distance dynamics. The former 
may be calculable, but the latter is completely nonperturbative. The strong phases for two-body 
nonleptonic B meson decays, derived in the QCDF or PQCD approaches, are all of the short-distance 
type. That is, the long-distance phases from final-state interaction (FSI) have been neglected in both 
approaches. It is then essential to look for any means to justify this assumption. In this subsection I 
review the recent study on this subtle and complicated topic. 

Most estimates of FSI effects in the literature |305] suffer ambiguities or difficulties. Some opinions 
favor small FSI effects in two-body nonleptonic B meson decays. For example, the smallness of FSI 
effects has been put forward based on the color-transparency argument |S7j . The renormalization-group 
analysis of soft gluon exchanges among initial- and final-state mesons j306j has indicated that FSI effects 
are not important in two-body B meson decays. For a limited number of decay channels, one can extract 
the FSI phases A5 directly from experiments. The phases are large for D meson decays |307j and small 
for B meson decays [308J: 

A5 = 80° ±7°; D — > K~tt + /T( 7i + , (196) 

f 11°; B -> D + n-/D Q n- 
A5 < I 16°; B -> D + p-/D°p- (197) 

[ 29°; B -> D* + 7r-/D*°7r-. 

implying that FSI phases diminish as the initial mass increases |309j . 

However, an opposite opinion was raised in P09 . Take the decay B + (bu) — ► K°(sd)7T + (du) as an 
example. The gluons exchanged between d and s (or d) are hard, but the gluons between the soft spec- 
tator u and s (or d) are not so hard. By simple kinematics, the invariant mass vhq ~ (Aqcd'^_b) 1/ ' 2 — 
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Figure 40: A model FSI phase. 



1.2 GeV for E u ~ Aq CD is in the middle of the resonance region of the sd channel. Then long-distance 
interactions cannot be ignored between K° and 7r + . 

Some attempts have been made to estimate FSI effects. Applying the time-reversal operation on 
the decay amplitude A n = (n out \0\B) , one obtains 

A n ^ (B\TOT~ l \n in ) . (198) 

Insert a complete set of out states, and employ the symmetry propery S nn i — S n i ni S nn i — (n ou ^\n! m ') . 
From the time reversal invariance of strong interaction, Eq. ()198|) becomes 

A n = Y, S nwK> , (199) 
n' 

for a T-even decay operator O. Subtracting the complex conjugate of A from both sides in Eq. (jl99J) 
and dividing it by i, the familiar relation emerges 

2\mA n = Y,tnniA* n , , (200) 

n' 

with t = (S — The above expression states that the strong phase is the sum of the contributions 
over all intermediate states n' . 

If approximating the intermediate states n' by two-body states, those connected to the final state 
n by Pomeron exchange dominate in the sum, and t nn > will be almost imaginary. The decay amplitude 
A n will be also almost imaginary, no matter what operator is responsible for the decay. Certainly, 
this approximation may not be justified. Without the two-body-state approximation, one can take 
advantage of the presence of many states. A statistical approach or a random approximation then helps 
pro] . Note that since A n / and S nn ' come from two different sources, weak and strong interactions, the 
phase of product S n i n A* n , for n' ^ n takes equally likely a positive or a negative value as n' is varied with 
n fixed. While A n / is related to A n (n' ^ n) by rescattering, there exist so many states that the influence 
of n on n! can be disregarded. Therefore, the phase of S nn iA* n , takes random values as n! varies. Under 
this approximation, a typical FSI phase in the B — > ttti and Ktt decays has been estimated to be about 
20° from the meson- meson scattering at 5 GeV |311j . The limitation of the random approximation is 
that one can only compute the statistically likely values of FSI phases as their standard deviations from 
zero, instead of predicting values of individual phases. 

FSI not only generates a strong phase but also changes a magnitude of amplitude. It has been 
argued that the decay amplitude could be enhanced by a FSI factor [312J, 

E = expf- I" 6 ^'\ ds') . (201) 
V 7r J Sth s' - m z B J 

In the approximation of two-body intermediate states, the FSI phase approaches ±90° at high energies. 
The enhancement or suppression effect is then very large due to a constructive integration over s' . In the 
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Figure 41: The enhancement (suppression) factor E vs E max . The lower curve (E < 1) for S(s') > 
shown here, and the upper curve for 5(s') — > —S(s'). 



random approximation, the sign of S(s') may fluctuate with s'. In this case, the effect of enhancement 
and suppression would be much smaller. An estimate for the B — > Kn, titt decays with one model FSI 
phase motivated by experiments has been performed. Choose the FSI phase as shown in Fig. 0U] 309J: 
S(s') rises to large values (~ 90°) around 2 GeV and falls linearly to zero at E max = 0{m B ). This S(s') 
suppresses the decay amplitude, as indicated by the lower curve in Fig. since the support of the 
phase integral is mostly below mg, where l/(s' — m 2 B ) is negative. A 10~40% correction to the decay 
amplitude is possible |3U9j . 

An application of Eq. (|199|) to the study of the B — > decays was provided in |293j . The matrix 
equation can be formally solved to give 

A = S 1/2 A° , (202) 

where A is an arbitrary real vector of the same dimension as A. One may consider the vector A as 
representing the decay amplitude in the absence of the final state phases from strong interaction. Assume 
elastic final state rescattering in the B — > Dtt modes. Assigning the real factorization amplitudes in 
FA to A , one obtains 



(203) 



where the superscript / denotes FA. 

To obtain S 1 / 2 , one derives the scattering matrix t through the optical theorem. This way makes 
transparent the mechanism involved in final-state rescattering [293 . By means of Eq. (|200|) and the 
various scattering mechanism displayed in Fig. H21 one has 



f A D o n - N 






\ 




| = S 1 ' 2 






\ A D o n o / 




l A f 


/ 



( r 0+K 



I 





r* + r* a 





r* + W*a + K) 



J 



(204) 



where r e , r a , and r$ parametrize the rescattering effects from charge exchange, annihilation, and flavor 
singlet (Pomeron) exchange, respectively, as defined in Fig. H21 The definition S — 1 + it then gives the 
constraints, 



r + r e 
r o + ^(3r a - r e ) 



2 sin 5 3/2 e i5 ^ 2 



2 sin Su 2 e 



/2 



(205) 



where 5i/ 2 and 5 3 / 2 are the phases associated with the corresponding isospin amplitudes. 

To perform the global fit, the data are put into the left-hand side of Eq. (J203|) . The amplitudes in 
FA with 02/^1 ~ 0.25 (input from the B — > J/tpK decay) and some model form factors are put into 
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Figure 42: Pictorial representation of (a) r e (charge exchange), (b) r a (annihilation) and (c) r$ (singlet 
exchange) . 



the right-hand side. The numerical analysis leads to [293 

(1 + iro)e- 2iS ^ 2 = 0.45 + 0.50i , 
ir ee -^/2 = 0.55 _ o.50i , 

ir a e~ 2iS3 ^ = 0.14 + 0.04i . (206) 

One then realizes the relative importance of the various rescattering mechanism: the Pomeron exchange 
and the charge exchange give roughly similar effect, larger than that from the annihilation mechanism. 
It was then concluded that FSI is crucial for enhancing a.2 from 0.25 for the B — > J/ipK decay to 0.5 
for the B° — > D 7t° decay p93j. Similar formalism has been applied to two-body charmless modes, and 
large FSI phases were also postulated. 

I emphasize that the major assumption in the above analysis is the absence of short-distance phases 
in two-body nonleptonic B meson decays. If short- distance phases exist, the decay amplitudes in Fig.l4*2l 
have carried phases already before rescattering. Therefore, the absence of short-distance phases must 
be assumed, so that the amplitudes before rescattering can be identified as the real amplitudes on the 
right-hand side of Eq. ()203|) . If short-distance phases do exist, the FSI effects could be small, and 
the conclusion on the relative importance of different mechanism might be altered. For example, the 
B — > Dtt branching ratios can be accounted for in the PQCD approach without resort to FSI as shown 
in Sec. O 

A stringent test of the small FSI assumption is provided by measuring the B — > KK decays. In 
particular, large observed B® — > K ± K T branching ratios and CP asymmetry in the B d — > K°K° modes 
will imply large FSI effects. So far, the PQCD predictions for the branching ratios [313 , 



B{B + - 


-> K + K°] 


= 1.47 x 10~ 6 


B(B~ - 


-> K-K°] 


= 1.84 x 10~ 6 


B(B° d - 


-> K+K~) 


= 3.27 x 10~ 8 


B(B° d - 


-> K-K + ) 


= 5.90 x 10~ 8 


B(B° d - 


-> K°K°) 


= 1.75 x 10~ 6 
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Figure 43: Intrinsic charm in the B-meson can mediate the decay to a strange, charmless final state via 
the weak transition b — > sec. 

B(B° d -> K°K°) = 1.75 x KT 6 . (207) 
are still below the experimental bounds. 

7.2 Intrinsic Charm and Charming Penguin 

One of the higher-power contributions comes from the Fock states of arbitrarily many particles. For 
example, a B meson bound state contains 

\B~) = ipbu\bu) + ip b ug\bug) + ip badd \budd) + ip bacS \bucc) H . (208) 

The Fock state decomposition is usually performed at equal light-cone time using light-cone quantization 
in the gauge A + = |314[ 1315] . The non- valence partons in the higher Fock states are generated by 
QCD splitting mechanism. All the partons in a Fock component are almost on mass shell with long 
lifetimes, and interact with each other through multiple infrared gluon exchanges. This is the reason 
they are intrinsic to the hadron structure. The intrinsic heavy quarks are part of the hadron bound 
state ]316j . Due to the hierachical structure of the CKM matrix elements, the weak transition b — > sec 
is doubly Cabibbo enhanced with respect to the b — > suu transition. This is the argument for the 
potential importance of the intrinsic charm (IC) |317j . In contrast, a perturbative correction to the weak 
transition matrix element can produce a cc pair through gluon splitting. The quark pair is generally not 
multiply connected to the partons of the bound state, and extrinsic to the hadron structure. Generally, 
"intrinsic" contributions in B meson decays are of higher twist, whereas "extrinsic" contributions are 
of higher order in a s . 

Some estimates based on phenomenological hints suggest that the IC probability in the B meson 
could be as large as a few percent |318j . The slight excess in the inclusive B — > J/ipX yield at low 
J/ip momentum |319j implies the presence of the B — > J/ipD^-*' channel, which occurs through IC, 
though such an effect could also be generated by the B — > J/ipKn decay 320J. IC could help to 
understand the large B —>■ rj'K , rj'X branching ratios. A valence cc component in the rj' meson has 
been introduced to resolve this discrepancy |321j . but the decay constant f^, ~ —2 MeV, defined via 

(0|c7 /1 75c|77 / (p)) = iffipv,, is too small |322j . Discussions on the above subjects can be found in |323j . 
It was proposed recently |324j that IC in the B meson could be probed by measuring the B —>■ J/i/jevX 
decay. 

As emphasized in |317j . IC plays an important role in the B — ► Kn decays. Considering only the 
valence contribution, the decay amplitude A(B° — > K + ir~) is written as 325J 

A(B° d - K+n-) = VusV^E, - P X GIM ) - VuVZPi . (209) 

The parameter E\ denotes the contribution from the W emission topologies, P\ denotes the penguin 
topologies, and ppiM con tains penguin contributions, which vanish in the m c = m u limit. Beyond the 
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Figure 44: Intrinsic charm in the B-meson can mediate the decay to a strange, charmless final state via 
the weak transition b — > sec. 

valence approximation, the additional contribution of IC through Fig. 03] arises, which is expressed as 
V cs V* b A 1 ^. There are no hard gluon exchanges across the weak vertex in Fig. 021 so that the computation 
of the hard scattering amplitude factorizes. One portion of the weak vertex mediates the annihilation 
of the be quarks, and the other describes the amplitude for the (cq) state to emerge with the parton 
content of the final state, namely sq'qq. Note that the above two pieces of the hard scattering 
kernels are in fact convoluted together over the momentum fraction of the light-cone wave function for 
the (bc)(cq) state, which is still unknown. 

The amplitude for the (cq) state to emerge as sq'qq cosmetically resembles Fig. !25Tf). Hence, the 
IC contribution has been parametrized as |317j 



where fs* ~ 0.317 GeV [326] arises from the annihilation of the be pair, and the remaining factors come 
from an estimate of the lower half of the diagram in Fig. 031 The factorizable annihilation amplitude 

has been calculated in the PQCD approach [22] . It is easy to get ai{mb/2) / a & {m b /2) ~ —20. The 
parameter B reflects the probability amplitude to find the B meson in an IC configuration, as well as 
an adjustment for the ~ 50% penguin enhancement [23]. Hence, one has B ~ 2(0.02)/3. The impact 
of IC on the B — * Ktt decays has been investigated in |317j based on Eq. (|210|) . and the results are 
exhibited in Fig. El Note that IC can act to either enhance or decrease the CP asymmetry. The IC 
contribution |y4* c |/|Pi| ~ O(10) % is a nontrivial fraction of the penguin parameter Pi, which reflects 
the accuracy we can reach in calculating the effective value of P\. 

"Charming penguins" refer to the nonperturbative piece of the extrinsic c quark loop, whose con- 
tributions are parametrized into a term P\. Theoretically, their importance is not clear. However, a 
recent analysis of the B — > Kir, tttt, Krj, and Krf decays indicates that they could be important. The 
parameter Pi includes not only the charming penguin contributions, but also annihilation and penguin 
contractions of penguin operators. It does not include leading emission amplitudes of penguin operators 
(O3-O6), which have been explicitly evaluated in FA. In this respect, it is a general parameterization 
of all the perturbative and nonperturbative 0(AQCT>/mb) contributions of the operators O5 and Oq. 
Pi has the same quantum numbers and physical effects as the original charming penguins proposed in 
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40.9±1.0 


3.70±0.42 


0.86±0.06±0.14 


230±30±15 


1.16±0.03±0.04 


Fk{MI) 
0.32 ±0.12 


18.6 ± 1.1 

J B(7T+7r-)xl0 6 


B(K + rr°) x 10 6 
11.5 ± 1.3 

B(tt + ii )x10 6 


B(K°tt + )x1Q 6 
17.9 ± 1.7 

B(n°TT°) x 10 6 


B(K°tt°)x1Q 6 
8.9 ±2.3 


0.27 ±0.08 


5.2 ±0.6 


4.9 ± 1.1 


<3.4 BaBar 





Table 10: Values of the input parameters used in the global fit. 



Mode 


UTA 

B (10- 6 ) \Acp\ 


\V u b/V C b\ 
B (10- 6 ) \A CP \ 


7T + 7r _ 
7T+7T 

ttV 


8.9 ±3.3 0.37 ±0.17 
5.4 ±2.1 
0.44 ±0.13 0.61 ±0.26 


8.7 ±3.6 0.39 ±0.20 
5.5 ±2.2 
0.69 ±0.27 0.45 ±0.27 


K+ti- 
K+n° 


18.4 ±1.0 0.21 ±0.10 
10.3 ±0.9 0.22 ±0.11 
19.3 ± 1.2 0.00 ±0.00 
8.7 ±0.8 0.04 ±0.02 


18.8 ±1.0 0.21 ±0.12 
10.7 ± 1.0 0.22 ±0.13 
18.1 ± 1.5 0.00 ±0.00 
8.2 ± 1.2 0.04 ±0.03 



Table 11: Predictions for CP-averaged branching ratios and absolute values of the CP asymmetries 
|Acp|. The left (right) columns show results obtained using constraints on the CKM parameters p and 
r] obtained from the unitarity triangle analysis (UTA) (the measurement of IV^/V^I). The last 
four channels are those used for fitting the charming penguin parameter P 1 . 



The charming penguin contribution in Pi can be, taking the B® — > K + tt~ decay as an example, 
included in the term Pi in Eq. (|209|) . Another 0(AqcdA^&) contribution from the GIM-penguin p^ IM 
has been neglected for simplicity |228j . The B — > Kir data do not constrain this parameter very 
effectively, since its contribution is doubly Cabibbo suppressed with respect to Pi. The remaining 
7r + 7r~ mode alone is not sufficient to fully determine the complex parameter P 1 G/A/ . However, the 
GIM-penguin contribution may be able to enhance the B(B — ► 7r°7r°) up to few xl0~ 6 |329j . 

Using the inputs collected in Table ITOl 228]. the complex parameter P\ = (0.13 ± 0.02) e ±J ( 114±35 )° 
in units of f. K F B ' K {m'^) has been determined from the data fitting, which has the expected size of 
0(AQCD/"ife)- Note that the sign of the phase is practically not constrained by the data. Table HT1 [228 
shows the predicted CP- averaged branching ratios and the absolute value of the CP-asymmetries \Acp\ 
for the B — > Kit and tttt modes. The angle 03 is determined through the effect of interefence terms 
in the B — > Kir branching ratios. The nonperturbative parameter Pi with an additional phase in the 
amplitudes makes the extraction of 3 difficult. It has been checked using the |\4b/K&|-constrained fit 
that almost any value of 03 is allowed, given the uncertainty on Pi. 

I stress that the intrinsic penguins and charming penguins represent distinct dynamics. However, 
both contributions are Cabibbo-enhanced and contain an 0(1) Wilson coefficient. Indeed, the parameter 
Af can also be absorbed into Pi and P^ 1M . Hence, a global fit can not distinguish them. Note that the 
intrinsic penguin contribution in |317j was parametrized to be proportional to the annihilation penguin 
amplitude calculated in the PQCD approach, such that its phase, including the sign, can be fixed. 

8 Conclusion 

In this article I have reviewed the two fundamental tools in QCD perturbation theory, collinear and 
kr factorization theorems, in which soft dynamics and hard dynamics of a process are factorized into 
hadron wave functions and hard kernels, respectively. Both factorization theorems can be constructed 
in a gauge-invariant way up to all orders, such that infrared-finite and gauge-invariant predictions 
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tractions becomes important, lo improve perturbative expansion, threshold resummation tor colhnear 
factorization and kx resummation for kr factorization can be performed. 

I then discussed the various QCD approaches to exclusive B meson decays available in the literature, 
which were developed based on the above fundamental concepts for perturbaton theory. I reviewed the 
recent progress made in these approaches, emphasizing the basic ideas behind and the comparison of 
their phenomenological implications. The competition among different approaches and the confronta- 
tion between theoretical predictions and experimental data have stimulated tremendous progress. In 
this section I briefly summarize the advantage of each method and the important issues, whcih require 
further investigation in order for a more solid and complete framework. 

The advantage of LCSR is that both soft and hard contributions can be analyzed in the same 
framework, and that it is easy to examine the self-consistency of expansions in a s and in \jm^. Recently, 
an essential step toward the evaluation of nonfactorizable contributions (from twist-3 hadron distribution 
amplitudes) has been made. However, the analysis is not yet complete. The urgent subject is to include 
the nonfactorizable contributions associated with the twist-2 distribution amplitudes, i.e., from the 
diagrams in Fig. These diagrams can not only moderate the scale dependence of predictions, but 
introduce strong phases, which are necessary for generating direct CP asymmetries. Since the strong 
phases come from higher-order contributions, compared to the leading soft transition form factors, they 
are expected to be small, the same as those from the QCDF approach. 

For QCDF, the advantage is its explicit factorization picture in the heavy quark limit. Talking 
about the treatment of leading contribtuions, it is most complete among all approaches, since the scale 
and scheme dependences have been greatly reduced. However, the end-point singularity appearing at 
subleading level (twist-3 nonfactorizable amplitudes and annihilation amplitudes) makes QCDF less 
predictive. This is also the reason the experimental constraint on the unitarity angles derived from 
QCDF is not very strong. Except for the end-point singularity, another challenge comes from the 
explanation of the possible large direct CP asymmetries in the two-body nonleptonic B meson decays. 

The PQCD approach, based on kx factorization theorem, is free of the end-point singularities. 
Therefore, the nonperturbative inputs are only universal hadron distribution amplitudes, without the 
transition form factors and the arbitrary infrared cutoffs due to the singularities. So far, the PQCD 
predictions are in agreement with data, and phenomenologically successful. However, it is still under 
debate whether the crucial Sudakov effect is strong enough to suppress the end-point contribution for 
the physcial mass mg ~ 5 GeV. The urgent subject is then to calculate higher-order corrections, and 
examine whether they converge. Before proving this, PQCD is not yet a self-consistent theory. The 
calculation can also verify the argument about the characteristic hard scale y/ Arris for exclusive B 
meson decays, which was made based on the hard-scattering picture. This characteristic scale is the 
key for PQCD to explain the B — > VP data. 

SCET provides the most systematic framework for constructing collinear factorization formulas of 
exclusive B meson decays at large recoil. Its advantage is that contributions characterized by different 
scales can be separated easily. A progress has been made recently in deriving the symmetry relations 
among various transition form factors, and the power corrections and radiative corrections to these 
relations. However, to make explicit predictions, which can be compared with data, it is necessary to 
calculate the Wilson coefficients of effective operators from the matching of the effective theory to the 
full one. The calculation is as complicated as in other approaches. 

The advantage of LFQCD is its simplicity. At leading power, the formula for a transition form 
factor involves only an overlap integral of the initial- and final-state hadron wave fuctions. Both soft 
dynamics and hard dynamics have been absorbed into the wave functions. The inclusion of nonvalence 
contributions has been worked out, which is crucial for guaranteeing the covariance of predictions. 
The next step is to extend the LFQCD to two-body nonleptonic decays, especially to nonfactorizable 
amplitudes. This extention will be a challenge. 

Though some study has been done, the behavior of the B meson bound state is still not clear. This 
subject is important, since the B meson wave functions, including those associated with higher Fock 
states (such as intrinsic charms), are the input of all QCD approaches. There is still no reliable method 
to analyze long-distance FSI effects in exclusive B meson decays. Even an ansatz does not exist, in 
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assumption, i.e., on the absence of short-distance phases, is not convincing as explained m Sec. 16.21 and 
in Sec. 17.11 Another crucial contribution from charming penguins also requires investigation. More 
clever ideas are necessary for understanding the above nonperturbative dynamics. 

As demonstrated in this article, exclusive B meson decays exhibit exciting QCD dynamics, which 
could be, fortunately, studied in a self-consistent way due to the large b quark mass. Some plausible 
mechanism has been explored. For example, annihilation contribution may not be as small as we 
thought with the hint from the possible large CP asymmetry in the B® — ► 7r + 7r _ decay. However, more 
challenging topics are waiting for our effort. 
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